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1  Introduction 


Sediment  Sampling 

All  studies  of  sedimentation  and  sediment  transport  require  some  degree  of  field 
data.  These  data  are  necessary  whether  the  study  consists  of  simple  observation  and 
correlation  of  empirical  evidence  or,  at  the  other  end  of  the  spectrum,  a  complex  three 
dimensional  numerical  model.  It  has  been  suggested  (USAGE,  1989)  that  a  staged 
approach  to  sediment  studies  is  the  best  approach  when  the  degree  of  the  problem  is 
unknown.  This  staged  approach  allows  ever  more  complex  steps  such  that  the  present 
step  builds  on  the  last.  Even  such  an  efficient  and  prudent  method  as  the  step  method 
requires  accurate  sediment  data. 

Necessary  data  to  conduct  sediment  studies  vary  as  to  the  complexity  of  the 
problem  and  the  approach,  but  for  alluvial  streams  can  generally  be  defined  as: 

*  Chaimel  geometry  (cross-section  size,  cross-section  shape,  reach  length,  water 
surface  or  energy  grade  line  slope,  extent  of  obstructions  and  width  of  flow) 

*  Stream  flow  characteristics  (i.e.,  discharge,  velocity,  depth  of  flow) 

*  Water  temperature 
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*  Suspended  sediment  characteristics  (grain  size  and  gradation  or  size  fraction 
distribution) 

*  Bed  sediment  characteristics  (grain  size  and  gradation  or  size  fraction 
distribution) 

It  can  quickly  be  reasoned  that  the  collection  and  recording  of  such  variables  would  be 
both  time  consuming  and  expensive.  The  question  now  becomes,  how  much  data  is 
necessary?  By  nature,  sediment  data  is  quite  variable.  Water  discharge  vs.  sediment 
discharge  relationships  are  typically  plotted  on  log-log  paper,  with  water  discharge  on  the 
abscissa  and  sediment  discharge  on  the  ordinate  axis.  For  a  given  water  discharge,  the 
associated  sediment  discharge  can  vary  over  several  log  cycles.  Furthermore,  other 
important  variables  such  as  grain  size  distribution  in  the  stream  bed  and  water  column  can 
vary  by  several  orders  of  magnitude.  Faced  with  these  degrees  of  variability,  the  historic 
answer  to  the  lead  question  has  been  to  collect  as  much  data  as  time  and  money  would 
allow  (USAGE  1948)  and  to  hope  that  the  data  covered  sufficiently  variable  conditions  to 
provide  representative  results  when  placed  into  the  transport  functions  and  used  as 
predictors.  This  approach  is  haphazard  at  best.  At  worst,  it  can  lead  to  erroneous 
conclusions  drawn  from  insufficient  or  misleading  data  sets. 
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Previous  Research 


Sediment  Data  Collection  and  Analysis 

Data  collection  in  general  tends  to  be  viewed  by  the  hydraulic  engineering 
community  as  somewhat  routine.  Perhaps  this  explains  why  little  is  written  in  the 
literature  on  this  subject  when  compared  to  other  subjects  which  are  undoubtedly  seen  as 
more  research  worthy.  However,  data  collection  in  general  is  always  a  significant,  if  not 
the  major  cost  for  hydraulic  and  sediment  studies.  Collection  and  analysis  of  sediment 
data,  in  particular,  tends  to  be  quite  expensive.  This  is  due  to  the  variability  of  the 
parameters  measured  (Guy,  1970).  The  majority  of  the  research  in  this  area  was 
conducted  by  a  fairly  limited  group  of  researchers  during  the  1940's  through  the  1970's. 
Research  literature  after  this  time  is  dominated  by  computer  applications  and  new 
numerical  methods  for  analyzing  and  predicting  sediment  movement. 

The  problem  of  data  collection  for  sediment  studies,  when  addressed  in  the 
literature,  typically  focuses  on  sediment  measurement  techniques.  These  sources,  of 
which  there  are  many,  relate  to  equipment  (Vanoni,  1975;  Shen,  1971;  Guy  and  Norman, 
1970;  Inter-Agency  Committee  on  Water  Resources,  Subcommittee  On  Sedimentation, 
Reports  Number  1  (1940),  2  (1940),  5  (1941),  6  (1952),  11(1957)  and  13  (1961))  and 
measurement  and  analysis  techniques  (Vanoni,  1975;  Shen,  1971;  Shen  and  Hikkawa, 
1980;  Inter-Agency  Committee  on  Water  Resources,  Subcommittee  On  Sedimentation, 
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Reports  Number  3  (1941),  4  (1941),  7  (1943),  8  (1948),  10  (1953),  and  12  (1957)). 

The  time  interval  for  taking  sediment  samples  has  been  discussed  (Guy  and 
Norman,  1970)  in  terms  of  adequacy  of  sediment  samples.  These  discussions  indicate 
that  adequate  sampling  may  require  that  samples  be  collected  on  daily,  hourly  or  even  bi- 
hourly  sampling.  However,  Guy  and  Norman  concede  that  each  sample  taken  costs  the 
same  at  a  station  and  that  monthly  sampling  over  a  long  period  may  yield  more 
information  and  be  more  cost  effective  than  daily  sampling  for  a  short  period.  The  fact 
that  daily  samples  can  often  yield  little  information  for  the  investment  has  been 
corroborated  (Martin,  1988).  It  has  further  been  suggested  that  samples  be  of  sufficient 
quantity  to  cover  the  natural  random  variation  and  si2e  gradation  of  transported  sediment 
(Porterfield,  1972). 

The  only  literature  identified  which  deals  with  the  number  of  samples  necessary  to 
construct  an  adequate  data  set  of  suspended  sediment  samples  is  the  work  of  H.  P.  Guy 
(1968).  The  technique  was  actually  intended  to  determine  if  sediment  samples  taken  in  a 
vertical  were  representative.  Its  purpose  was  to  replace  the  necessity  for  computing  the 
algebraic  means  of  different  populations  in  order  to  quantify  a  population.  Guy 
developed  a  statistical  technique  based  on  the  deviation  of  measurements  fi-om  the  mean 
of  the  suspended  sediment  concentration.  The  procedure  would  assess  the  quality  of  a 
mean  for  a  group  of  data  or  the  quality  of  a  coefficient  based  on  a  ratio  of  sample  means. 
He  presented  this  technique  as  a  nomograph  which  would  either  evaluate  an  existing  data 
set  or  predict  the  number  of  samples  required  based  on  an  assumed  sum  of  squared 
deviations  from  the  mean  and  a  desired  level  of  confidence. 
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In  summary,  previous  research  as  to  how  many  sediment  samples  are  necessary  to 
achieve  the  best  results  has  been  very  limited.  What  research  that  has  been  done  has 
almost  exclusively  focussed  on  equipment,  measurement  and  analysis  techniques.  Guy 
(1968)  developed  some  guidance  for  the  number  of  samples  require  to  provide  a  desired 
level  of  confidence  in  the  population  mean  of  samples  in  a  vertical  section  of  a  steam 
cross  section.  Since  Guy’s  work,  there  has  been  no  additional  research  as  to  how  many 
sediment  samples  are  necessary  to  describe  the  sediment  characteristics  at  a  sediment 
sampling  station.  It  is  currently  assumed  that  sediment  samples  will  be  collected  for  a 
given  duration,  with  this  duration  controlled  primarily  by  the  funds  available  for  the 
study. 

Sediment  Transport  Functions 

Sediment  transport  functions,  per  se,  (the  term  ‘functions’  is  used  throughout 
because  as  a  group  these  exist  in  the  form  of  graphical  solutions,  complex  formulations 
and  regression  derived  equations)  are  not  the  topic  of  the  research  presented  herein. 
However,  since  they  are  used  as  part  of  a  procedure  for  determining  the  adequacy  of 
sediment  data  at  a  sampling  station,  it  seemed  prudent  to  provide  some  background  as  to 
their  theoretical  development.  Sediment  transport  in  streams  occurs  in  three  ways.  It  is 
either  moved  along  and  in  contact  with  the  bed,  by  saltation  (skipping  and  bouncing 
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along  the  bed)  or  in  suspension.  One  or  more  of  these  may  occur  simultaneously.  The 
transport  processes  of  rolling  and  saltation  occur  veiy  near  the  bed.  Sediment  transported 
in  these  manners  is  denoted  bed  load  (Vanoni,  1975).  The  zone  of  bed  load  transport  is 
variously  reported  as  two  bed  particle  diameters  thick  (Mien,  1995)  or  a  few  bed  particle 
diameters  (USAGE,  1989).  The  suspended  load  is  sediment  transport  that  occurs  in  the 
zone  above  the  zone  of  bedload  transport  and  extends  to  the  water  surface.  Total  load  is 
the  combination  of  bed  load  and  suspended  load.  To  complete  this  definition  of  terms, 
the  suspended  load  is  further  divided  into  bed  material  load  and  wash  load.  Bed  material 
load  is  comprised  of  sediment  particles  found  in  appreciable  quantities  in  the  bed.  Wash 
load  is  comprised  of  very  small  particles  found  in  very  small  quantities  in  shifting 
portions  of  the  bed  (Vanoni,  1975). 

Sediment  sampling,  as  conducted  by  the  United  States  Geological  Survey  in  the 
St.  Francis  River  basin,  collects  samples  of  suspended  sediment  from  the  entire  water 
column  with  the  exception  of  the  portion  within  0.3  to  0.4  feet  of  the  bottom  (Guy  and 
Norman,  1970).  This  unsampled  zone  contains  all  of  the  bed  load  in  sand  bed  streams 
and  a  small  potion  of  the  suspended  load. 

The  pioneering  work  in  sediment  transport  computation  was  done  by  M.  P.  Duboy 
(1879).  His  work  was  based  on  the  concept  that  sediment  moves  in  thin  layers  along  the 
stream  bed.  The  basis  of  his  approach  was  that  the  applied  bed  shear  stress,  Tq,  must 
exceed  the  critical  shear  stress,  t^,  in  order  to  initiate  motion.  Duboy’s  equation  is 
presented  below,  in  English  units. 


6 


(1) 


^bv 


-0.0125 -0.019rfJ 
d~^ 


where 

qbv  =  volume  of  bed  material  in  motion,  per  foot  of  width 
ds  =  particle  size  in  millimeters 
To  =  boundary  shear  stress  in  Ib/ft 

Mien  (1995)  has  classified  sediment  transport  functions  into  three  categories 
based  on  the  basic  approach  involved  in  their  derivation.  These  are  (1)  formulations 
based  on  advection-diffusion,  (2)  formulations  based  on  energy  concepts,  and 
(3)  empirical  methods  based  on  regression  analysis  and  graphical  procedures.  While 
many  sediment  transport  functions  have  been  proposed  since  the  pioneering  work  of 
Duboys,  no  one  function  accurately  and  consistently  works  for  all  applications.  In 
practice,  one  typically  selects  several  functions  based  on  their  applicability  to  a  given 
stream.  This  applicability  is  based  on  the  similarity  of  data  from  which  the  function  is 
developed  to  that  found  on  the  stream  in  question.  The  function  results  are  then 
compared  in  some  manner  to  field  data  and  the  one  that  best  fits  the  field  data  is  selected 
(Mien,  1995). 

Three  approaches  are  discussed  below.  These  represent  one  from  each  type  of 
formulation,  as  presented  above.  Each  of  these  was  used  in  this  dissertation. 
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Colby’s  Method.  Colby  (1964)  used  the  very  complex  Einstein  bed  material  function 
(Einstein,  1950)  which  built  on  Einstein’s  earlier  bedload  transport  function  (Einstein, 
1942).  Einstein’s  function  was  based  on  the  statistical  probability  of  a  particle  initiating 
motion.  He  combined  this  with  advection-dififusion  principles  to  develop  his  bed  material 
function.  Colby  developed  his  relationship  based  on  an  immense  amoimt  of  stream  and 
flume  data.  He  used  measurements  from  a  score  of  streams  including  the  Niobrara  River, 
the  Middle  Loup  River,  the  Colorado  River  and  the  Mississippi  River.  In  addition,  he 
augmented  these  data  with  flume  data  collected  in  three  different  size  flumes.  Colby 
applied  the  Einstein  bed  material  function  to  the  computation  of  bed  material  discharge 
for  a  given  size  of  bed  sediment  over  a  wide  range  of  velocities,  depths  and  water 
temperatures  to  fill  gaps  in  the  stream  and  flume  data.  He  also  added  a  correction  for 
sediment  concentrations  in  excess  of  10,000  parts  per  million  and  for  temperatures  more 
or  less  than  60  degrees  Fahrenheit.  Colby  presented  his  results  as  a  series  of  four  graphs. 

In  order  to  apply  the  Colby  method,  one  need  only  know  water  depth,  mean  velocity, 
temperature  and  sediment  concentration.  Colby’s  method  is  theoretically  applicable  over 
a  range  of  water  depths  from .  1  feet  to  1 00  feet  deep,  velocities  from  1 .0  to  1 0.0  feet  per 
second,  median  bed  sand  sizes  from  0.10  to  1.0  millimeters  in  diameter,  water 
temperatures  from  32-100  degrees  Fahrenheit  and  sediment  concentrations  up  to 
200,000  parts  per  million.  Its  wide  range  of  applicability,  ease  of  use  and  relative 
accuracy  have  all  contributed  to  the  popularity  of  this  method.  It  was  particularly  popular 
prior  to  the  advent  of  computer  based  sediment  transport  calculations.  In  recent  years,  it 
has  lost  some  prestige  to  other  methods  that  provide  better  results  and  whose  complexity 
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are  not  an  obstacle  to  computer  generated  results  (Y ang,  1991). 

Yang’s  Method.  Yang  (1973)  seized  upon  the  findings  of  Cook  (1945)  and  Bagnold 
(1960)  to  develop  his  method  which  is  founded  on  the  notion  of  stream  power.  Stream 
power,  as  was  proposed  by  Cook  and  Bagnold,  is  the  rate  at  which  a  stream  loses  energy 
per  unit  area  of  the  boundary.  For  three-dimensional  flow  and  a  boundary  that  includes 
the  banks,  the  stream  power  is  given  as 

yRSV 

or  the  product  of  the  total  shear  and  the  mean  velocity.  Yang  modified  this  definition  and 
expressed  stream  power  as  rate  of  potential  energy  dissipated  per  unit  weight  of  water. 
This  is  defined  as  VS  or  the  product  of  the  mean  velocity  and  the  energy  slope.  Yang’s 
dimensionless  unit  stream  power  equation  is: 


d  u 

logC^^^  =5.435  -0.286IogO)— -0.4571og-  ^ 


0) 


d  u 

+[  1.799  “0 .4091oga)—  -0.3 1 41og— ] 

V  O) 

vs 

.  Iog[— — ^] 

(0  (0 


(2) 


where 

Cpp„  =  the  concentration  of  bed-material  discharge  in  parts  per  million  by 
weight. 
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0)  —  the  average  fall  velocity,  in  feet  per  second  of  sediment  particles  of 
diameter  djo 

ds  =  the  particle  size  in  feet  at  which  50  percent  of  the  bed  material  by  weight 
is  finer 

u  =  the  inematic  viscosity  in  feet  squared  per  second 
u*  =  the  shear  velocity  in  feet  per  second 
V  =  the  average  velocity  in  feet  per  second 
S  =  the  energy  slope  in  feet  per  foot 
Vc  =  the  average  flow  velocity  in  feet  per  second  at  incipient  motion 
The  sediment  transport,  in  tons  per  day,  is  related  to  Yang’s  concentration  value  by 

Qs  =  Qw  C  k  (3) 


where 

Qs  =  sediment  transport  in  tons  per  day 

Qw  =  water  discharge  in  cubic  feet  per  second  (English  units)  or  cubic  meters 
per  second  (metric) 

C  =  sediment  concentration  in  parts  per  million 
k  =  0.0027  (English)  or  0.0864  (metric) 

The  coefficients  for  Yang’s  equation  were  developed  from  463  sets  of  laboratory  data 
Yang’s  equation  was  developed  from  data  ranging  from  0.137  to  1.71  millimeters  median 
sieve  diameter  and  for  water  depths  of  0.037  to  49.9  feet.  Yang’s  equation  has  proven  to 
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be  remarkably  reliable  over  a  wide  range  of  sand  bed  streams  in  actual  use.  In  a 
comparison  with  seven  other  widely  accepted  transport  equations,  based  on  1,1 19  sets  of 
laboratory  data  and  319  sets  of  stream  data,  Yang’s  equation  proved  most  accurate  (Yang 
and  Wan,  1991).  However,  no  one  equation  has  yet  been  developed  that  is  universally 
more  accurate  than  all  others  (Mien,  1995). 

Brownlie’s  Method.  Brownlie  (1981)  chose  yet  a  different  approach  from  the  two 
previously  discussed.  He  reasoned  that  sediment  data  from  streams  was  too  variable  to 
lend  itself  to  physics-based  analysis.  He  therefore  developed  his  method  based  on 
dimensional  analysis  and  regression  techniques.  Brownlie  proposed  that  a  relationship 
for  sediment  concentration  should  have  the  general  form 


C=/[?,5,g,p,v,p^,c/5o,oJ 


where 


C  =  sediment  concentration 
q  =  water  discharge,  per  unit  width 
S  =  water  surface  slope 
p  =  density  of  water 
u  =  kinematic  viscosity 
Ps  =  specific  gravity  of  the  sediment  particle 
djo  =  median  grain  size 
Og  =  gradation  coefficient 
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Brownlie  rearranged  the  eight  independent  variables  into  five  dimensionless  groups, 
shown  below. 


C=Aq, 


P,-P, 


(5) 


where 

q.  =  dimensionless  water  discharge 
R  =  Reynold’s  number 

and  the  other  variables  are  as  described  above. 

By  calculation  and  further  combining  of  the  terms,  Brownlie  arrived  at  the  final  form  of 
his  equation. 


r  P  “P 

c=y[F  ^,5.0  .i? 


(6) 


where 

Fg  =  grain  Froude  number 

r  =  hydraulic  radius 

Rg  =  grain  Reynolds  number 

and  the  other  variables  are  as  described  above. 
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Using  a  data  base  of  unprecedented  size  (5,263  laboratory  sets  and  1,764  stream  sets), 
Brownlie  used  multiple  regression  analysis  to  develop  his  relationship.  The  final  form  is 


C=7115cJF  -F  )i  ”8^  0.6601  Jl_ 

g  go'  , 

^50 


-0.3301 


(6) 


where 

Cp  =  coefficient  which  equals  1.0  for  laboratory  data  and  1.268  for  stream  data 

Fgo  =  criticzd  grain  Froude  number 

and  the  other  variables  are  as  described  above. 

Brownlie’s  method  is  applicable  to  rivers  with  sand  size  bed  material  ranging  from 
0.02-76.113  millimeters  median  diameter,  average  velocities  between  0.472  -  5.188  feet 
per  second,  depths  of  .06  -  53.9  feet  deep,  slopes  of  .0000021  -  .0126  feet  per  foot  and 
water  temperatures  of  32  -  96.8  degrees  Fahrenheit. 

Entropy  and  Water  Resources. 

The  notion  of  using  entropy  as  a  measure  of  the  degree  of  ignorance  as  to  the  true 
state  of  a  thermodynamic  system  dates  to  the  19th  century  (Boltzman,  1 872).  It  was  then 
that  Boltzman  described  entropy  as: 


H=k  log  p 


(7) 
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where  H  is  the  entropy,  p  is  the  probability  of  system  state,  and  k  is  Boltzman’s  constant. 
C.  E.  Shannon  incorporated  entropy  into  his  Mathematical  Theory  of  Communications 
(Shaimon,  1948a,  1948b).  This  work  was  developed  in  five  parts  and  included 
23  Theorems  associated  with  the  transmission,  receipt  and  recovery  of  messages  sent 
electronically.  In  the  first  of  these  parts,  “Discrete  Noiseless  Systems”  Shannon 
discussed  communication  theory  as  relates  to  discrete  systems  such  as  telegraphs.  He 
introduced  entropy  as  a  measure  of  uncertainty  involved  between  a  signal  sent  and  the 
one  received.  He  stated  that  the  ratio  of  the  entropy  of  a  source  to  the  maximum  entropy 
it  could  have  should  be  denoted  its  relative  entropy.  This  value  was  identified  as  the 
maximum  compression  possible  when  signals  are  encoded  for  a  given  language’s 
alphabet.  One  minus  the  relative  entropy  was  defined  as  the  redundancy.  The 
redundancy  of  ordinary  English  was  stated  to  be  50%.  This  meant  that  written  English  is 
one  half  determined  by  the  structure  of  the  language  and  one  half  is  freely  chosen  by  the 
author.  Part  II  addressed  a  discrete  channel  with  noise  introduced  which  corrupted  the 
sent  signal.  Shannon  developed  mathematical  theories  to  determine  errors  on  the 
receiving  end  and  correct  these  errors  based  on  a  corrective  chaimel  capacity  greater  than 
the  average  error,  determined  partially  through  entropy  calculations.  Part  III  outlined  the 
mathematical  preliminaries  of  the  next  parts.  In  this  part,  the  notion  of  discrete 
probability  distributions  was  introduced.  Discrete  distributions  represent  continuous 
functions  as  groups  of  occurrences,  such  as  those  associated  with  the  dots  and  dashes  of 
telegraphy.  Continuous  distributions,  such  as  those  associated  with  the  spoken  word, 
were  also  presented.  Of  note  was  the  finding  that,  in  the  discrete  case,  entropy  measures 
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in  an  absolute  way  randomness  of  the  chance  variable.  In  the  continuous  case,  the 
entropy  is  a  measure  of  randomness  relative  to  an  assumed  standard,  namely  the 
coordinate  system.  However,  in  either  case  entropy  is  an  important  measure  of 
uncertainty  due  to  the  fact  that  the  derived  concepts  of  information  rate  and  channel 
capacity  depend  on  the  difference  of  two  entropies  and  this  difference  does  not  depend  on 
the  coordinate  system  as  each  of  these  two  term’s  coordinate  systems  are  changed  by  the 
same  amount.  The  final  parts  addressed  the  capacity  of  a  continuous  channel  and  fidelity 
evaluation  factors. 

More  current  research  using  entropy  concepts  developed  by  Shannon  has  been 
presented  by  E.  T.  Jaynes,  (1957a,  1957b)  in  the  field  of  informational  theory  and 
statistical  mechanics  for  the  prediction  of  equilibrium  thermodynamic  properties.  In 
these  articles,  Jaynes  introduced  informational  theory  which  provides  a  constructive 
criterion  for  setting  up  probability  distributions  on  the  basis  of  partial  knowledge.  This  is 
developed  to  a  type  of  statistical  inference  which  is  called  the  maximum  entropy  estimate 
which  is  formally  known  as  the  Principle  of  Maximum  Entropy.  This  estimate  provides 
the  least  biased  estimate  possible  based  on  the  information  given.  Statistical  mechanics 
were  presented  as  a  form  of  statistical  inference  rather  than  a  physical  theory.  Jaynes 
considered  the  terms  “entropy”  and  “uncertainty”  as  synonymous.  He  sketched 
Shannon’s  proof  that  the  quantity  which  is  positive,  increases  with  increasing  uncertainty 
and  is  additive  for  independent  sources  of  uncertainty  is  given  by: 
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where 


(8) 


HiPy..p„)  =  -kY^  plnp. 

»=i 


H  is  entropy 

K  is  a  positive  constant 

Pi  is  the  probability  distribution 

The  above  equation  is  the  Shannon  Entropy  Function.  Jaynes  sketched  the  proof  of  this 
function  as  follows.  “The  variable  x  can  assume  the  discrete  values  (x, , ..  jcJ.  Our  partial 
understanding  of  the  processes  which  determine  the  value  of  x  can  be  represented  by 
assigning  corresponding  probabilities  (p^  We  ask,  with  Shannon,  (1948a)  whether 

it  is  possible  to  find  any  quantity  H(p,.../7j  which  measures  in  a  unique  way  the  amoimt 
of  uncertainty  represented  by  this  probability  distribution.  It  might  at  first  seem  very 
difficult  to  specify  conditions  for  such  a  measure  which  would  ensure  both  uniqueness 
and  consistency,  to  say  nothing  of  usefulness.  Accordingly,  it  is  a  very  remarkable  fact 
that  the  most  elementary  conditions  of  consistency,  amounting  to  really  to  only  one 
compositional  law,  already  determines  the  function  H(p,...;7j  to  within  a  constant  factor. 
The  three  conditions  are: 

(1)  H  is  a  continuous  function  of  the/?; . 

(2)  If  all/?;  are  equal,  the  quantity  A{ri)  =  H(l/n,...,l/«)  is  a  monotonic  increasing 
function  of  n. 

(3)  The  composition  law.  Instead  of  giving  the  probabilities  of  events  (x,  ...x„) 
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directly,  we  might  group  the  first  k  of  them  together  as  a  single  event,  and  give  its 
probability  w,  =  (p,  +.../?*!;  then  the  next  m  possibilities  are  assigned  the  total  probability 
W2  =  (p*+i+—+A+<«)5  etc.  When  this  much  has  been  specified,  the  amount  of  uncertainty  as 
to  the  composite  events  is  H(w,...w,.).  Then  we  give  the  conditional  probabilities 
(p,/w„...,/7/wi)  of  the  ultimate  events  (X|..jc*),  given  that  the  first  composite  event  had 
occurred,  the  conditional  probabilities  for  the  second  composite  event,  and  so  on.  We 
arrive  ultimately  at  the  same  state  of  knowledge  as  if  the  had  been  given  directly, 

therefore  if  our  information  measure  is  to  be  consistent,  we  must  obtain  the  same  ultimate 
uncertainty  no  matter  how  the  choices  were  broken  down  in  this  way.  Thus,  we  must 
have 


^(P\-Pn)  =  H(w,...W2)  +  w,  H  (p,/w„...,Vwi) 

+  W2H(p*^,/W2,...,Pk+„/W2)  +... .  (9) 

The  weighting  factor  w,  appears  in  the  second  term  because  the  additional  uncertainty 
H  (p,/w„...,/7/w,)  is  encountered  only  with  probability  w,.  For  example, 

HC/2,  1/3, 1/6)  =  H(‘/2,  1/2)  +  1/2H(  2/3,  1/3). 

From  condition  (1),  it  is  sufficient  to  determine  H  for  all  rational  values 

Pr «/!«,.  (10) 

with  Kj  integers.  But  then  condition  (3)  implies  that  H  is  determined  already  from  the 
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symmetrical  quantities  A(n).  For  we  can  regard  a  choice  of  one  of  the  alternatives 
as  a  first  step  in  the  choice  of  one  of 

n 

/=i 

equally  likely  alternatives,  the  second  step  of  which  is  also  a  choice  between  k,  equally 
likely  alternatives.  As  an  example,  with  n  =  3,  we  might  choose  («„«2>«3)  =  (3,4,2).  For 
this  case  the  composition  law  becomes 

In  general,  it  could  be  written 

+Y,pA(n)  =A(^p,)  (12) 


In  particular,  we  could  choose  all  «,•  equal  to  m,  where-upon  (12)  reduces  to 

A(m)+A(n)=A(mn)  (13) 

Evidently  this  is  solved  by  setting 

A(n)=K  Inn  (14) 

where,  by  condition  (2),  K>0.  For  a  proof  that  (14)  is  the  only  solution  of  (13),  we  refer 
the  reader  to  Shannon’s  paper  (1948a).  Substituting  (14)  into  (12),  we  have  the  desired 
result.” 
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H(pj...pJ=K  ln(ln)-KYp,lnn. 

=  -Kl,p\np. 


(15) 


The  application  of  entropy  principles  to  water  resources  has  been  applied  by 
Sonuga  (1972, 1976)  concerning  the  rainfall-runoff  relationship  and  hydrologic  frequency 
analysis.  Amorocho  and  Espildora  (1973)  used  the  Principle  of  Maximum  Entropy  to 
analyze  the  uncertainty  in  stream  flow  computed  using  the  Stanford  Watershed  Model. 

V.  P.  Singh  has  been  one  of  the  most  prolific  researchers  to  apply  entropy  principles  to 
water  resources  in  recent  years.  His  works  have  addressed  the  derivation  of  frequency 
distributions  commonly  used  for  hydrologic  analyses  using  Jayne’s  Principle  of 
Maximum  Entropy,  (Singh,  Rajagopal,  and  Singh,  1986).  A  stochastic  model  for  total 
sediment  yield  predictions  based  on  direct  runoff  volume  using  entropy  principles  was 
developed  and  compared  favorably  for  several  watersheds,  (Singh  and  Krstanovic,  1987). 
V.  P.  Singh  and  M.  Fiorentino  (1992)  published  an  exhaustive  historical  perspective  of 
entropy  applications  in  water  resources.  The  period  covered  included  1962-1991 . 

Subjects  covered  included  fluvial  geomorphology  (Leopold  and  Langbein,  1962);  open 
channel  hydraulics  (Chiu,  1986, 1988;  Barbe  et  al.,  1991);  redimdancy  measures  for  water 
distribution  networks  (Awumah  et  al.,  1 99 1 );  assessment  of  numerical  model 
performance  (Amorocho  and  Espildora,  1973;  Chapman,  1986);  data  acquisition  related 
to  sampling  intervals  for  ammonia  concentrations  in  France  and  stream  flow 
measurements  in  Turkey  (Harmancioglu,  1981, 1984);  rainfall  network  design 
(Krstanovic  and  Singh,  1988a,  1988b);  catchment  modeling  (Jowitt,  1991);  and  a  large 
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body  of  work  addressing  the  use  of  entropy-based  probability  distribution  function 
parameter  estimation  and  comparisons  of  this  method  with  other  parameter  estimating 
methods  (Jowitt,  1979;  Singh  and  Jain,  1985;  Jain  and  Singh,  1986;  Arora  and  Singh, 
1987a,  1987b,  1989;  Singh,  Singh,  1985a,  1985b,  1985c,  1988;  Arora  and  Singh,  1989; 
Fiorentino  et  al.,  1987a,  1987b;  Singh  et  al.,  1989, 1990a,  1990b;  and  Singh  and 
Rajagopal,  1987). 


Knowledge  Gap 


There  exists  a  paradox  in  sedimentation  engineering  of  the  latter  20th  Century. 
Resources  to  collect  data  and  conduct  sediment  studies  are  limited.  However,  the 
repercussions  of  improper  or  ill  advised  engineering  have  never  been  greater,  due  to 
increased  land  costs,  dense  population  and  the  litigious  nature  of  society. 

The  advent  of  high  speed  computing  and  software  to  take  advantage  of  this 
capability,  have  led  to  unprecedented  ability  to  analyze  and  predict  sediment  behavior 
under  a  wide  variety  of  conditions.  The  area  where  knowledge  is  lacking  is  in  die  data 
collection  arena.  As  stated  previously,  these  data  are  expensive  to  collect  and  analyze. 
Field  data  are  necessary  to  calibrate  and  verify  complex  numerical  models.  Data 
collected  in  good  faith  can  still  lead  to  erroneous  analytical  conclusions  due  to  the 
variable  nature  of  the  data. 
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What  is  needed  is  a  definitive  answer  to  the  question,  “how  much  data  are 
necessary?”,  and  its  corollary  “Do  1  have  enough  data?”.  If  one  knew  how  to  determine 
the  answer  to  these  questions,  then  the  level  of  confidence  in  the  analytical  solutions 
could  be  determined.  Sediment  sampling  programs  could  then  be  designed  to  address 
these  concerns  in  the  long  view.  Existing  data  sets  could  be  used  with  a  known  degree  of 
confidence  in  the  short  view. 

Report  Outline  and  Organization 

The  major  objective  of  this  research  is  to  offer  specific  guidance  as  to  how  much 
sediment  data  are  needed  to  provide  the  maximum  amount  of  information  for  use  in 
describing  the  sediment  characteristics  of  a  fluvial  system.  This  is  accomplished  by 
introducing  a  data  set  upon  which  the  research  is  based,  describing  the  development  of  a 
procedure  to  meet  the  objective,  and  executing  said  procedure.  Findings  and  results  are 
then  presented  and  conclusions  drawn  as  to  the  success  in  meeting  the  objective. 

The  report  is  organized  in  six  chapters  as  follows.  The  source  and  history  of  the 
data  set  used  in  this  research  are  presented  in  Chapter  2,  along  with  a  discussion  as  to  its 
quality  and  applicability  to  this  research.  Data  base  construction,  data  manipulation,  and 
initial  computations  necessary  to  the  research  are  presented  in  Chapter  3.  Statistical 
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analyses  and  entropy  calculations  are  presented  in  Chapter  4.  Development  and 
description  of  a  procedure  for  evaluating  sediment  data  sets  as  to  their  adequacy  for 
selecting  a  representative  sedinient  transport  function  and  determining  if  they  are  of 
optimum  length  is  presented  in  Chapter  5.  Chapter  6  contains  conclusions  as  to  the 
usefulness  of  the  procedure,  examples  of  potential  applications  for  the  procedure  and 
recommendations  for  additional  research. 
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2  Research  Data 


Background 


General 

The  St.  Francis  River  Basin,  located  in  southeast  Missouri  and  eastern  Arkansas, 
in  the  early  1900's  was  a  vast  expanse  of  bottom  land  hardwoods  and  poorly  drained, 
swampy  areas.  In  1906,  a  report  was  published  (USAGE,  1906)  that  outlined  a  plan  for 
drainage  of  this  region  and  conversion  to  cropland.  In  general,  the  plan  and  its  later 
revisions  called  for  construction  of  a  flood  control  reservoir  in  the  upper  basin.  The 
remainder  of  the  plan  called  for  massive  channeli2ation  of  the  tributaries  of  the 
St.  Francis  River,  construction  of  a  levied  floodway  along  the  course  of  the  river,  and 
realignment  of  the  river  itself.  The  present  day  St.  Francis  River  basin  drains 
5,159  square  miles.  The  basin  is  about  215  miles  long  and  has  a  maximum  width  of 
53  miles.  The  basin  slopes  vary  from  3  feet  per  mile  in  the  upper  basin  to  0.5  feet  per 
mile  in  the  lower  basin.  The  bed  material  is  largely  sand  with  an  average  grain  size  of 
1.0  millimeter  in  diameter.  As  called  for  in  the  various  plans,  the  basin  has  been  highly 
developed  hydraulically  with  structural  features  including  a  dam,  siphons,  inverted 


23 


siphons,  levied  flood  ways,  channel  cutoffs,  extensive  drainage  of  wetlands  along  its 
tributaries  and  improved  channel  reaches  (USAGE  1985).  The  reservoir,  Wappapello 
Dam  and  Lake,  controls  drainage  from  1,3 10  square  miles,  or  about  25  percent  of  the 
basin  area. 


Sediment  Sampling  Program 


Changing  channel  conditions  and  land  use  during  the  period  1930-1975  created  a 
highly  dynamic  river  system  that  was  in  a  constant  state  of  flux  and  thus  a  system-wide 
equilibrixim  was  not  possible.  As  the  system  sought  to  reach  an  new  equilibrium 
condition,  numerous  sediment-related  problems  arose.  These  included  but  were  not 
limited  to  scour  holes  at  bridges,  bank  failures,  unraveling  of  tributary  channels, 
excessive  channel  sedimentation,  and  overbank  deposition  of  sediments.  In  order  to 
assist  in  operating  and  maintaining  this  complex  network  of  hydraulic  works,  the 
United  States  Army  Corps  of  Engineers,  Memphis  District  decided  to  enter  into  an 
agreement  with  the  U.S.  Geological  Survey  (USGS)  to  collect  sediment  and  water  quality 
data  at  various  sites  within  the  basin. 

In  1977,  the  U.S.  Army  Engineer  District,  Memphis  Hydraulics  Branch  proposed 
a  sediment  sampling  program  to  create  a  data  base  of  sufficient  depth  and  breadth  to 
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address  these  myriad  sediment  problems.  The  author  assisted  the  Hydraulics  Branch 
Chief  in  designing  the  initial  sediment  sampling  program  and  maintained  the  records  and 
performed  various  analyses  using  the  data  during  the  period  1977-1983. 

The  objective  of  the  sampling  program  was  to  collect  sufficient  data  to  define  the 
sediment  budget  for  the  basin,  while  providing  data  necessary  for  engineering  analyses 
aimed  at  addressing  sediment  problems  within  the  basin.  Twenty-four  sites  were 
originally  identified  as  necessary  to  define  the  sediment  budget  for  the  basin.  Fourteen  of 
these  were  on  the  main  stem  of  the  St.  Francis  River  or  its  man-made  floodway  and  the 
remainder  on  major  tributaries.  These  were  arranged  such  that  important  tributary 
sediment  sources  could  be  identified.  In  addition,  sediment  transport  by  reach  could  be 
computed  fi'om  the  measurements,  depicting  areas  of  likely  scour  or  deposition.  The  data 
collected  went  far  beyond  the  basic  sediment  data  needed  to  define  a  crude  sediment 
budget.  Both  bed  material  and  suspended  sediment  samples  were  collected  and  analyzed. 
These  data  provided  the  necessary  detailed  input  to  support  accurate  regime  and 
sedimentation  studies  using  the  data  collected.  The  collection  of  data  at  the  various  sites 
was  to  be  done,  on  a  cost  reimbursable  basis,  by  the  U.S.  Geological  Survey  office  in 
Little  Rock,  AR.  The  cost  for  data  collection  at  each  site  averaged  $5,000  for  a  year  or 
about  $420  for  each  measurement.  The  program  was  initiated  in  October,  1977  and  data 
were  collected  on  a  water-year  basis,  that  is,  from  October  through  the  following 
September.  The  original  24  site  locations  are  shown  in  Figure  1 . 
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POPLAR  BLUFF 


ST.  FRANCIS  RIVER 
NEAR  GLENNONVILLE 

WILHEMINACUT-OFF 


ST.  FRANCIS  RIVER 
ATRSK 

ST.  FRANCIS  RIVER 
NEAR  POWE 


ST.  FRANCIS  RIVER 
AT  ST  FRANCIS 


LOCUST  CREEK 
DITCH  (HIGHWATER) 


BIG  SLOUGH 
DITCH 


LOCUST  CREEK  DITCH 


THOMPSON  CREEK 


JONESBORO 

ST.  FRANCIS  RIVER 
AT  LAKE  CITY 


ST.  FRANCIS  RIVER 
FLOODWAY  DITCH  60 


STRAIGHTSLOUGH 

WYNNE 

CLARK  CORNER  CUT-OFF 

FORREST  CITY 


ST.  FRANCIS  RIVER 
NEAR  PIGGOTT 

KENNETT 


COCKLE  BURR  SLOUGH 


RIGHT  HAND  CHUTE  OF 
UTTLE  RIVER  AT  RIVERVALE 


MEMPHIS 


ST.  FRANCIS  RIVER 
AT  MADISON 


Figure  1.  Original  24  sediment  sampling  stations,  1977 


Previous  uses  of  these  data  include  but  are  not  limited  to  providing  boundary 
conditions  for  sediment  analyses;  looking  at  comparative  transport  rates  for  different 
reaches  of  the  stream  system;  to  investigate  sediment  contributions  from  the  various 
tributaries;  and  to  address  claims  of  sand  deposition  on  farm  lands. 


Sampling  Method 


The  Equal  Transit  Rate  (ETR)  method  was  used  to  collect  the  samples.  In  this 
method,  samples  are  obtained  at  equally  spaced  verticals  across  the  stream.  The  sampler 
is  lowered  and  raised  at  a  constant  rate.  Since  the  vertical  velocity  distribution  is  not 
uniform,  the  sample  obtained  in  this  manner  is  quantitatively  weighted  according  to  the 
velocity  through  which  it  passes.  Such  a  sample  is  considered  to  be  discharge  weighted 
because  with  a  uniform  transit  rate  each  increment  of  velocity  or  discharge  in  the  vertical 
is  given  the  same  amount  of  time  to  enter  the  sampler.  This  equal  spacing  between 
verticals  and  equal  transit  rate,  both  up  and  down  in  all  verticals,  yields  a  gross  sample 
proportional  to  the  total  stream  flow  (Guy  and  Norman  1970).  The  number  of  verticals 
sampled  varied  from  station  to  station  and  even  at  a  given  station  with  varying  discharge; 
however,  the  number  of  verticals  was  typically  near  twenty.  This  is  noted  as  adequate  for 
most  ail  streams  (Guy  and  Norman,  1970).  In  this  method,  the  transit  rate  for  the 
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sampling  bottle  is  set  by  the  vertical  with  the  highest  discharge.  The  round  trip  is  timed 
so  that  the  sample  bottle  is  just  filled  at  the  end  of  the  upward  trip.  The  samples  were 
composited  to  produce  a  single  suspended  sediment  concentration  and  material  for  a 
suspended  sediment  sample  grain  size  distribution  curve.  Likewise,  bed  material  samples 
were  collected  at  the  verticals  and  composited  to  produce  one  sample  for  which  a  bed 
material  grain  size  distribution  curve  was  developed. 

Discharge  measurements  were  made  in  conjunction  with  the  collection  of  the 
sediment  samples.  These  measurements  were  primarily  made  by  the  USGS.  At  seven  of 
the  24  stations,  the  Corps  of  Engineers  were  responsible  for  the  discharge  measurements 
and  continue  to  do  so.  The  USGS  timed  their  sampling  to  coincide  with  the  discharge 
meeisurement  dates  at  these  seven  locations. 


Data  Collected 


At  each  station  the  following  37  articles  of  data  were  noted,  collected  or 
computed; 

*  Measurement  number 

*  Individual  making  computations 

*  Individual  checking  computations 
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*  Station  number 


*  Date 

*  Party  members 

*  Width  of  stream 

*  Area  of  flow 

*  Average  velocity 

*  Gage  height 

*  Discharge 

*  Method  of  measuring  velocity  (either  0.6  depth  or  0.2  and  0.8  averaged) 

*  Number  of  verticals 

*  Gage  change  during  measurement 

*  Meter  data  (T5qpe,  number  coefficients  etc.) 

*  Method  of  sampling  (  wading,  boat,  etc.) 

*  Gage  readings 

*  Rating  for  measurement 

*  Flow  conditions 

*  Weather 

*  Air  temperature 

*  Water  temperature 

*  Flow  control 

*  Whether  or  not  overbank  flow  was  present 


*  Remarks  by  party 


*  Width  of  vertical 


*  Depth  at  vertical 

*  Velocity  observation  depth 

*  Meter  revolutions 

*  Time  for  revolutions 

*  Velocity  at  measurement  point 

*  Mean  velocity  in  vertical 

*  Adjustments  for  horizontal  angle  of  meter  line 

*  Area  in  vertical 

*  Discharge  in  vertical 

*  Suspended  sediment  gradation 

*  Bed  material  gradation 


Laboratory  Data  Analysis 


The  data  collected  were  sent  to  a  USGS  laboratory  for  analysis.  These  analyses 
yielded: 

*  Suspended  sediment  concentration,  in  milligrams  per  liter  or 
parts  per  million  (PPM) 
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*  Suspended  sediment  discharge  in  tons  per  day 

*  Suspended  sediment  diameters  (using  the  grade  scale  proposed 
by  the  American  Geophysical  Union  (Lane,  1947)) 

*  Bed  sediment  diameters  (Lane,  1947) 

The  percent  of  the  sample  contained  in  each  size  class  was  determined  using  sieve 
screens  for  the  early  data  and  by  use  of  a  Visual  Accumulation  Tube  (VAT)  shortly  after 
sampling  began.  These  were  the  methods  from  which  grain  size  distributions  were 
determined  for  each  bed  and  suspended  sediment  sample.  The  sieve  data  particle  sizes 
were  referred  to  as  sieve  diameter  whereas  the  VAT  data  particle  sizes  were  referred  to  as 
fall  diameter.  These  data  are  presented  in  Appendix  A  under  the  sample  S  AS  output. 

The  following  code  was  used  since  the  software  used  would  not  allow  a  numerical 
designation  of  a  variable. 

Suspended  sediment  size,  determined  on  the  basis  of  seive  diameter. 

SSSA  =  2.0  mm 
SSSB  =  1.0  mm 
SSSC  =  0.5mm 
SSSD  =  0.25  mm 
SSSE  =  0.125  mm 
SSSF  =  0.0625  mm 

Suspended  sediment  size,  determined  on  the  basis  of  fall  diameter. 

SSFA  =  2.0  mm 
SSFB=  1.0  mm 
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SSFC  =  0.5  mm 


SSFD  =  0.25  mm 
SSFE  =  0.125  mm 
SSFF  =  0.0625  mm 

Bed  sediment  size,  detennined  on  the  basis  of  seive  diameter. 


BMSA  = 

32.0  mm 

BMSB  = 

16.0  mm 

BMSC  = 

8.0  mm 

BMSD  = 

4.0  mm 

BMSE  = 

2.0  mm 

BMSF  = 

1.0  mm 

BMSG  = 

0.5  mm 

BMSH  = 

0.25  mm 

BMSI  = 

0.125  mm 

BMSJ  = 

0.0625  mm 

Bed  sediment  size,  determined  on  the  basis  of  fall  diameter. 
BMFA  =  2.0  mm 
BMFB=  1.0  mm 
BMFC=  0.5  mm 
BMFD  =  0.25  mm 
BMFE  =  0.125  mm 


BMFF  =  0.0625  mm 


Data  Selection 


The  data  set  collected  for  use  in  this  research  consisted  of  all  sediment  data 
collected  in  the  St.  Francis  River  Basin  for  Water  years  1978-1988  inclusive.  This 
extensive  data  set  was  analyzed  to  determine  how  much  of  the  data  set  could  or  should  be 
included  in  the  analysis.  Criteria  were  established  as  follows; 

(1)  The  stations  selected  should  be  representative  of  different  flow  regimes  so  that 
relationships  developed  would  not  be  station  specific.  Ideally,  there  should  be  at  least 
one  station  from  the  upper  basin,  one  from  the  middle  basin,  and  one  from  the  lower 
basin. 

(2)  The  flow  should  be  largely  confined  to  the  main  channel.  Sediment  samples 
were  composited  at  the  stations  and  one  analysis  performed.  Limiting  the  stations  to  ones 
that  had  main  channel  flow  only  would  remove  potential  bias  introduced  by  a  large 
number  of  overbank  samples  being  added  to  the  composite  sediment  samples. 

(3)  Tributary  stations  were  not  considered.  These  records  were  suspect  due  to 
influences  of  backwater  from  the  St.  Francis  and  often  very  low  or  no  flow  was  recorded 
for  these  stations. 

(4)  It  must  be  possible  to  compute  any  needed  variables,  notably  the  slope,  for 
any  station  selected. 

(5)  The  record  should  extend  for  the  entire  eleven  years. 
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The  original  twenty-four  data  collection  stations'  identification  numbers,  names, 
frequency  and  chronology  of  sampling,  and  latitude  -  longitude  coordinates  are  given  in 
Table  1.  It  can  be  seen  from  Table  1  that  beginning  with  water  year  1983,  the  data 
collection  program  began  a  series  of  modifications.  In  this  year,  four  stations  located  on 
tributaries  were  deleted  due  to  poor  location  or  the  fact  that  their  contributing  basins  were 
of  such  a  size  that  the  monthly  collection  schedule  typically  recorded  "no  flow".  Fifteen 
stations  were  modified  such  that  data  were  collected  only  during  the  months  of  November 
through  Jime.  Only  five  of  the  original  stations  were  continued  as  originally  conceived 
with  monthly  collection  of  data.  These  decisions  were  economically  driven  as  the  cost  of 
the  sampling  program  was  under  scrutiny  by  upper  management.  Again,  beginning  with 
water  year  1985,  the  data  collection  program  was  further  modified  for  economic 
considerations.  Five  more  stations  were  deleted,  four  on  tributaries  and  one  on  the  main 
stem  of  the  St.  Francis  River. 

Based  on  the  above  selection  criteria,  a  station-by-station  analysis  was  conducted. 
Criteria  (3)  and  (5)  eliminated  1 1  stations  (stations  4,8,9,1 1,12,13,14,15,16,17,23). 

Criteria  (2)  eliminated  an  additional  five  stations  (6,7,10,18,19).  This  left  eight  possible 
stations  which,  unfortunately,  were  located  in  two  clusters,  one  in  the  upper-middle 
portion  of  the  basin  (stations  24,22,21,20)  and  one  in  the  lower  portion  of  the  basin 
(stations  1,2,3,5).  Two  stations  were  selected  from  the  upper  grouping.  Fisk  was  a  good 
selection  as  the  uppermost  station  in  the  basin.  It  had  a  good,  stable  cross  section,  little 
overbank  flow  and  the  necessary  slope  calculations  could  be  made  using  the  Glennonville 
station.  Fisk  was  thus  selected  for  analysis.  Stations  20,21  and  22  are  located  on  a  reach 
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Table  1 

Data  Station  Identification,  Chronology  and  Location,  St.  Francis  River  Basin, 

MO.  &  AR. 

Water  Year/  #  Samples  per  Water  Year 

CE  ID 

USGS  ID 

78  83  85 

Latitude/Longitude 

1 

07040000  St.  Fr.  R. 

35  46’  50" 

@  Madison,  AR. 

12  12  12 

90  12’  80" 

07047904  Clark  Corner 

35  08’  41" 

Cutoff  nr.  Colt,  AR. 

12  8  8 

90  39’  23" 

3 

07047900  St.  Fr.  Bay 

35  15’  34" 

@  Riverfront,  AR. 

12  12  12 

90  40’  48" 

4 

07047882  Straight  Slough 

35  21’  45" 

nr.  Birdeye,  AR. 

12  8  — 

90  39’  26" 

5 

07047815  Cross  County 

35  21’  38" 

Ditch  nr.  Birdeye,  AR. 

12  8  8 

90  39’  00" 

6&7 

07047810  St.  FR.  R. 

35  32’  15" 

nr.  Marked  Tree  AR. 

12  8  8 

90  29'  05" 

8 

07046600  Rt.  Hand  Chute 

35  40"  20" 

Little  R.  @  Rivervale,  AR. 

12  12  12 

90  29*  12" 

9 

07040496  Cockle  Burr 

35  51’  39" 

Slough  nr.  Monette,  AR. 

12  8  8 

90  19’  49" 

07040450  St.  Fr.  R. 

35  49'  16" 

@  Lake  City,  AR. 

12  12  12 

90  25*  56" 

11 

07040440  Thompson  CR. 

35  52’  11" 

nr.  Lester,  AR. 

12  —  — 

90  25'  56" 

12 

07040445  Big  Bay  Ditch 

35  52'  12" 

nr.  Lester  AR. 

12  —  — 

90  27'  45" 

13 

07040428  Locust  Cr.  Ditch 

35  55’  12" 

nr.  Dixie,  AR. 

12  —  — 

90  25’  56" 

14 

07040424  Eight  Mile  Creek 

35  59’  16" 

nr.  Paragould,  AR. 

12  8  — 

90  25'  39" 

(Continued) 
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Table  1  (Concluded) 


Water  Year/  #  Samples  per  Water  Year 


CE  ID 

USGS  ID 

78 

83 

85 

Latitude/Longitude 

15 

07040424  Locust  Cr.  Ditch 

35  58’  06” 

nr.  Paragould.  AR. 

12 

8 

— 

90  24'  17" 

16 

07040350  Big  Slough  Ditch 

36  02’  25” 

nr.  Paragould.  AR. 

12 

8 

— 

90  21'  39" 

17 

07040150  Varney  R. 

36  08’  20" 

nr.  Senath,  MO. 

12 

— 

90  13'  54" 

18 

07040130  St.  Fr.  R. 

36  14'  11" 

@  Holly  Island,  AR. 

12 

8 

8 

90  07’  52" 

19 

07040110  St  Fr.R. 

36  23' 50" 

nr.  Piggott,  AR. 

12 

8 

8 

90  04' 40" 

20 

070400100  St.  Fr.  R. 

36  27’  21’’ 

@  St.  Francis,  AR. 

12 

8 

12 

90  08'  13" 

21 

07040070  Wilhelmina 

36  30’  53" 

Cutoff  nr.  Campbell,  MO. 

12 

8 

8 

90  09*  30" 

22 

07040060  St.  Fr.  R. 

36  39'  38" 

nr.  Glennonville,  MO. 

12 

8 

8 

90  11’ 06" 

23 

07040057  St.  Fr.  R. 

36  39’  38" 

nr.  Powe,  MO. 

12 

8 

— 

90  08’  32" 

24 

07040000  St.  Fr.  R. 

36  46’  50" 

@  Fisk,  MO. 

12 

12 

12 

90  12’  08" 
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of  the  river  approximately  10  miles  in  length.  Due  to  their  proximity,  one  was  deemed 
representative  of  this  reach  of  the  river.  Of  the  three,  station  21,  Wilhelmina  Cutoff,  is 
the  most  interesting  from  a  sediment  transport  aspect.  It  is  located  on  a  man-made  reach 
that  reduced  the  river  length  by  some  7  miles.  It  experiences  very  high  rates  of  sediment 
transport  and  thus  would  provide  a  severe  test  for  any  method  developed  in  the  basin. 
Slope  calculation  could  be  made  using  station  22  observations.  Therefore,  Wilhelmina 
was  also  selected  for  inclusion  in  the  analysis. 

Of  the  lower  basin  stations,  station  1,  Madison,  is  located  in  the  backwater 
influence  of  the  Mississippi  River.  It  was  therefore  removed  from  consideration  for 
selection  or  use  to  compute  slopes  at  other  stations.  The  remaining  three  stations  were  all 
equally  likely  candidates  located  on  a  fourteen  mile  long  reach  of  the  river.  However,  of 
the  three,  station  2  offered  the  largest  drainage  area,  since  it  was  below  the  large  Straight 
Slough  tributary,  while  the  other  two  were  located  above  this  tributary.  Therefore, 
station  2,  Clark  Comer  Cutoff,  was  chosen.  The  three  stations  that  were  left  after 
consideration  of  the  criteria  offered  three  different  cases  with  which  to  conduct  research. 
They  represented  a  range  of  drainage  areas  from  1,370  square  miles  for  Fisk  to  over 
5,000  square  miles  at  Clark  Comer.  Peak  observed  water  and  sediment  discharges  range 
from  30,800  cubic  feet  per  second  and  36,341  tons  per  day,  respectively,  at  Clark  Comer 
to  10,400  cubic  feet  per  second  and  1,954  tons  per  day  at  Fisk. 

In  summary,  stations  2,  (Clark  Comer  Cutoff  near  Colt,  Arkansas),  21 
(Wilhelmina  Cutoff  near  Campbell,  Missouri),  and  24  (Fisk,  Missouri)  were  selected  for 
detailed  analysis.  These  stations  offered  varied  hydraulic  and  sediment  conditions  and 
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very  high  quality,  professionally  collected  data,  with  which  to  perform  the  research.  The 
station  locations  are  shown  in  Figure  2. 
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Figure  2.  Location  of  sediment  sampling  stations  used  in  the  analyses 


Research  Data  Base  Creation 


The  data  selected  for  analysis  were  gathered  from  several  different  sources.  Some 
of  the  data  collected  were  published  in  standard  USGS  format  (USGS  1977-1988).  A 
sample  of  these  are  shown  for  the  Fisk,  Missouri,  station  in  Figure  3.  While  usefiil,  the 
data  in  this  form  would  require  cumbersome  transfer  to  a  data  base  to  allow  numerical 
manipulation.  In  discussions  with  the  USGS  in  Little  Rock,  Arkansas,  it  was  discovered 
that  the  data  were  available  in  digital  format.  This  format,  referred  to  as  'Flat  File 
Format’  by  the  USGS,  contained  ail  the  data  shown  in  Figure  3. 

Additional  data  were  available  on  the  field  data  sheets.  These  sheets  are  forms 
filled  out  by  the  survey  party  members  on  site  while  actually  collecting  the  The 
information  recorded  is  that  described  in  the  'Data  Collected'  section  above.  A  copy  of  a 
sample  field  data  sheet  for  the  Fisk,  Missouri,  station  is  shown  in  Figure  4.  From  these 
sheets,  values  for  the  top  width  of  flow,  area  of  flow,  average  flow  velocity  and 
maximum  depth  of  flow  were  obtained.  In  addition,  average  depth  was  computed  by 
dividing  the  flow  area  by  the  top  width.  The  digital  data  set  was  thus  augmented  by 
adding  these  additional  values  by  means  of  hand  editing  the  files.  Additional  data  were 
computed  as  will  be  discussed  in  Chapter  3. 
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ST.  FRANCIS  RIVER  BASIN 


07040000  ST,  FKANCIS  RIVER  AT  FISK.  HO 
PERIOD  OF  RECORD.-- October  1977  Co  current  year. 


Bucler-Stoddard  County  Itnc 


WATER  QUALITY  DATA,  WATER  YEAR  OCTOBER  1982  TO  SEPTEMBER  1983 


DATE 

TIME 

MEDIUM 

AGENCY 

COL¬ 

LECTING 

SAMPLE 

(CODE 

HUMBER) 

(00027) 

AGENCY 

ANA¬ 

LYZING 

SAMPLE 

(CODE 

HUMBER) 

(00028) 

STREAM- 

FLOW, 

INSTAN¬ 

TANEOUS 

(CFS) 

(00061) 

SPE¬ 

CIFIC 

OON- 

DUCT- 

•AKCE 

(UMHOS) 

(00095) 

PH 

(STAND¬ 

ARD 

UNITS) 

(00400) 

TEMPER¬ 
ATURE 
(DEC  C) 
(00010) 

TRANS¬ 

PAR¬ 

ENCY 

(SECCHI 

DISK) 

(M) 

(00078) 

OCT 

19... 

MOV 

1230 

0 

80513 

80513 

335 

174 

7.7 

17.5 

.70 

15... 

29.-. 

JAM 

1400 

1150 

0 

0 

80513 

80513 

80513 

80513 

328 

1020 

172 

182 

7.9 

7.6 

9.0 

8.5 

.50 

.80 

10... 

FEB 

1400 

0 

80513 

80513 

10400 

85 

7.0 

8. -5 

.40 

14... 

MAR 

1250 

0 

80513 

80513 

3360 

203 

7.8 

6.5 

.90 

15... 

APR 

1430 

9 

80513 

80513 

1320 

204 

8.5 

12.0 

.50 

13... 

MAY 

1630 

9 

80513 

80513 

2360 

176 

8.0 

12.0 

.80 

09... 

JUN 

1645 

9 

80513 

60513 

7000 

108 

7.2 

17.5 

.50 

14... 

JUL 

0830 

9 

80513 

80513 

910 

160 

7.4 

22.0 

.40 

12... 

AUG 

1300 

9 

80513 

80513 

178 

179 

7.7 

29.0 

.20 

08... 

SEP 

1230 

9 

80513 

80513 

no 

223 

8.2 

28.5 

.20 

13... 

1400 

9 

80513 

80513 

61 

229 

8.5 

27.0 

.20 

DATE 

TIME 

OXYGEN, 

DIS¬ 

SOLVED 

(MG/L) 

(00300) 

OXYGEN, 

DIS¬ 

SOLVED 

(PER¬ 

CENT 

SATUR¬ 

ATION) 

(00301) 

BARO¬ 

METRIC 

PRES¬ 

SURE 

(MM 

OF 

‘MG) 

(00025) 

SEDI¬ 

MENT. 

SUS¬ 

PENDED 

(MG/L) 

(80154) 

SEDI¬ 

MENT, 

DIS¬ 

CHARGE. 

SUS¬ 

PENDED 

(T/DAY) 

(80155) 

SED. 

SUSP. 

FALL 

DIAM. 

X  FINER 
THAN 
.062  MM 
(70342) 

SED. 

SUSP. 

FALL 

DIAM. 

X  FINER 
THAN 
.125  MM 
(70343) 

SED. 

SUSP. 

FALL 

DIAM. 

X  FINER 
THAN 
.250  MM 
(70344) 

OCT 

19... 

HOV 

1230 

7.8 

82 

760 

33 

30 

95 

98 

100 

15.. . 

29.. . 
JAN 

1400 

1150 

11.6 

12.6 

101 

108 

760 

760 

31 

56 

27 

154 

66 

66 

75 

97 

94 

99 

10... 

FEB 

1400 

11.2 

97 

750 

38 

1070 

80 

90 

95 

14... 

MAR 

1250 

12.6 

103 

756 

49 

445 

68 

76 

97 

15... 

APR 

1430 

10.7 

101 

750 

29 

103 

97 

97 

98 

13... 

MAY 

1630 

9.3 

87 

760 

85 

542 

76 

90 

99 

09... 

JUN 

1645 

9.6 

101 

760 

31 

586 

78 

'84 

95 

14... 

JUL 

0830 

6,5 

75 

758 

71 

174 

94 

97 

100 

12... 

AUC  ■ 

1300 

6  *4 

84 

760 

58 

26 

88 

69 

90 

08... 

SEP 

1230 

7.5 

97 

757 

29 

8.6 

92 

96 

100 

13... 

1400 

7.9 

100 

760 

26 

4.3 

97 

98 

98 

Figure  3.  Sample  of  the  published  USGS  data  used  in  the  analyses  (Continued) 


ST.  FSAHCIS  RIVER  BASIH 

07040000  ST.  FRAKCZS  RIVER  AT  FISR,  HOo-OOKTlNUEO 
WATER  QUALITY  DATA.  WATER  YEAR  OCTOBER  1982  TO  SEPTEMBER  1983 


SED. 

BED 

BED 

BED 

BED 

BED 

'  BED 

SUSP. 

MAT. 

MAT. 

MAT. 

MAT. 

MAT. 

HAT. 

FALL 

FALL 

FALL 

FALL 

FALL 

FALL 

FALL 

OIAK. 

OZAH. 

DIAH. 

DZAK. 

DIAM. 

DIAH. 

DIAH. 

X  FINER 

Z  FINER 

X  FINER 

X  FINER 

X  FINER 

Z  FINER 

X  FINER 

DATE 

THAN 

THAN 

THAN 

THAN 

TUAN> 

THAN 

THAN 

TIKE 

.500  HK 

.062  KK 

.125  KK 

.250  KM 

•500  KM 

1.00  KM 

2.00  KM 

(70345) 

(80158) 

(60159) 

(60160) 

(80161) 

(80162) 

(80163) 

OCT 

19... 

1230 

— 

2 

3 

77 

99 

too 

_ 

KOV 

15... 

1400 

100 

2 

2 

79 

99 

100 

29... 

JAM 

1150 

100 

3 

4 

56 

98 

100 

— 

10... 

FEB 

1400 

100 

16 

27 

86 

100 

— 

- 

14... 

1250 

100 

2 

7 

83 

too 

KAR 

15... 

APR 

1430 

100 

3 

3 

38 

97 

100 

— 

13... 

HAY 

1630 

100 

2 

2 

55 

98 

100 

— 

09... 
JUN  « 

1645 

100 

8 

8 

15 

83 

100 

14...  . 

JUL 

0830 

— 

4 

6 

55 

90 

100 

- 

12... 

1300 

100 

4 

18 

96 

99 

100 

AUC 

08... 

SEP 

1230 

— 

3 

19 

95 

99 

100 

- 

13... 

1400 

100 

5 

18 

94 

99 

99 

100 

Figures.  (Concluded) 
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Figure  4.  Sample  field  data  sheet 


Data  Base  Operations 


In  order  to  provide  ease  of  manipulation  and  the  capability  to  perform  various 
numerical  functions  on  the  data,  it  was  deemed  necessary  to  use  a  software  program  or 
programs  that  could  manipulate  the  data,  perform  numerical  functions  on  the  data  and 
present  the  results  in  tabular  and/or  graphical  form.  In  the  course  of  the  analysis,  several 
programs  were  ultimately  used  to  accomplish  the  construction  of  and  operations  with  the 
data  base.  Initially,  all  of  the  data  were  input  into  the  Statistical  Analysis  Software  (SAS) 
Package  developed  by  the  SAS  institute,  Cary,  North  Carolina  (SAS,  1988).  A  sample  of 
the  output  from  this  software  is  shown  in  Appendix  A.  The  data  printed  in  this  format 
was  easy  to  read  and  reference  to  dates  or  other  values.  While  the  software  contained 
powerful  statistical  routines,  the  software  required  renewal  on  an  annual  basis  for 
continued  use.  For  this  reason,  other  software  was  used  for  statistical  analyses  as 
discussed  below. 

SuperCalc4  spreadsheet  software  (Computer  Associates,  1987)  developed  by 
Computer  Associates,  San  Jose,  California,  was  used  to  compute  stream  slopes  for  each 
observation  used  in  the  analyses.  The  gage  heights  were  input  for  adjacent  gages.  These 
were  added  to  the  gage  zero  and  the  difference  divided  by  the  distance  between  the  gages 
to  compute  the  water  surface  slope  for  each  observation.  The  slope  could  not  always  be 
calculated  for  each  raw  observation  due  to  missing  gage  data  at  either  of  the  two  stations 
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used  in  the  computations.  Therefore,  the  total  number  of  observations  was  reduced  for 


this  reason  by  twenty-five  at  Fisk,  two  at  Wilhelmina  and  five  at  Clark  Comer.  An 
example  of  this  output  is  given  in  Appendix  A. 

The  U.S.  Army  Corps  of  Engineers  Hydraulic  Design  Package  for  Channels 
(SAM)  (Thomas  et  al.,  1992)  was  used  to  compute  the  sediment  transport  for  the 
observations.  These  data  then  became  the  basis  for  additional  research  as  is  described  in 
chapter  3.  Grain  size  information  on  the  bed  material  is  required  by  the  SAM  program. 
The  sediment  samples  were  either  lost  or  not  taken  for  a  number  of  the  observations, 
rendering  sediment  transport  computations  impossible.  For  this  reason,  the  munber  of 
data  observations  available  for  analysis  was  reduced  by  two  at  Fisk,  two  at  Wilhelmina 
and  sixteen  at  Clark  Comer.  An  example  of  the  SAM  input  and  output  are  included  in  . 
Appendix  A.  Table  2  below  summarizes  the  data  adjustment  made  and  described  above. 


Table  2 

Data  Reduction 

Station 

Raw  Observations 

Reduction-Missing 

Reduction-  Missing 

Adjusted  Number  of 

Slope 

Bed  Samples 

Observations 

Fisk 

128 

25 

2 

101 

Wilhelmina 

103 

0 

2 

101 

Clark 

105 

5 

16 

84 

STATGRAPHICS  software,  (Manugistics,  1993)  developed  by  Manugistics, 


Rockville,  Maryland  was  used  for  statistical  analyses,  distribution  fitting  and  plotting  of 


45 


results.  The  input  files  for  this  software  were  extracted  from  those  used  in  the  S  AS 
software,  with  observation  deleted  where  no  slope  or  sediment  data  were  available,  as 
discussed  above.  This  software  was  used  to  perform  the  regression  and  distribution 
fitting  exercises  presented  in  Chapter  4.  Many  of  the  figures  in  this  report  were 
developed  directly  from  this  software.  A  sample  input  file  is  shown  in  Appendix  A. 

Quattro  Pro  V  6.0,  (Wordperfect,  1993)  developed  by  the  Novell  Applications 
Group,  WordPerfect,  Orem,  Utah,  was  used  for  the  entropy  calculations  presented  in 
Chapter  4  and  plotting  of  certain  results.  The  necessary  input  was  extracted  from  the 
other  data  sources  and  input  in  spreadsheet  format.  Some  of  the  figures  in  this  report 
were  developed  directly  from  this  software.  A  sample  input  file  is  shown  in  Appendix  A. 
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3  Initial  Approach  to  the  Problem 


Introduction 


The  first  step,  after  consolidating  the  data  base,  was  to  check  the  goodness  of  the 
samples  in  the  available  data.  Furthermore,  it  was  decided  that  the  data  for  the  Fisk, 
Missouri  station  would  be  used  to  develop  the  initial  methodology.  This  station  was 
selected  because  it  had  the  narrowest  range  of  water  flows  and  sediment  discharge. 
Subsequent  development  would  proceed  in  a  downstream  direction,  with  ever  higher 
water  flows  and  increasing  sediment  discharge.  None  of  the  data  were  arbitrarily 
eliminated  as  being  non-representative,  as  it  was  all  representative  of  data  sets  that  might 
reasonably  be  expected  to  be  analyzed  in  the  future.  Twenty-seven  of  the  total 
observations  at  the  Fisk  Station  were  eliminated  for  lack  of  slope  or  bed  material  data  as 
summarized  in  Table  2.  It  should  be  noted  that  the  basin  characteristics  are  sufficiently 
different  from  the  vicinity  of  the  Fisk  station  to  those  at  the  lower  end  as  to  provide  a 
fairly  rigorous  test  of  any  procedure  developed. 
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Sediment  Transport  Computations 


Computations  of  sediment  transport  were  necessary  for  each  observation.  These 
would  be  the  basis  for  determining  goodness  of  the  data  collected  and  to  assist  in 
determining  the  number  of  samples  required  to  adequately  describe  the  sediment 
transport  regime. 

Initially,  the  sediment  transport  computations  were  made  for  each  observation 
using  the  CORP  system  (Jones,  1977).  However,  this  method  required  input  for  each 
individual  observation  and  for  each  transport  function  investigated.  In  addition  to 
gathering  the  individual  hydraulic  input  parameters,  such  as  average  velocity,  kinematic 
viscosity,  depth  and  width  of  flow,  volume  of  flow,  etc.,  detailed  sediment  data  were  also 
required.  These  were  in  the  form  of  either  the  djo  or  the  d^j  bed  material  grain  size.  To 
determine  these  values,  it  was  necessary  to  plot  the  bed  material  data  as  a  grain  size 
distribution  curve.  An  example  of  this  is  shown  in  Figure  5.  This  was  clearly  too  slow 
and  cumbersome  for  the  large  number  of  solutions  needed.  To  address  this  problem,  the 
SAM  system  (Thomas  et  al.,  1992)  mentioned  above  was  used.  The  SAM  system  allows 
simultaneous  computation  of  sediment  transport  using  thirteen  sediment  transport 
functions.  These  include  functions  developed  by  Toffaleti  (Toffaleti,  1969);  Yang  (Yang, 
1973);  Ackers-White  (Ackers  and  White,  1973);  Colby  (Colby,  1964);  Toffaleti- 
Schoklitsch  (in  preparation);  Meyer-Peter  Muller  (Meyer-Peter  and  Muller,  1948); 
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Brownlie  (Brownlie,  1981);  Toffaleti-Meyer-Peter  Muller  (in  preparation);  Laursen- 
Madden  (Madden,  1985);  Laursen-Copeland  (Copeland  and  Thomas,  1989).  There  were 
also  three  modified  transport  functions.  Those  that  used  only  the  single  djo  grain  size  in 
the  computations  are  denoted  “djo.”  These  include  Yang,  Ackers-White  and  the  Meyer- 
Peter  Muller  equations.  The  other  versions  of  these  transport  functions  have  been 
modified  by  The  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  Hydraulics 
Laboratory  (Thomas  et  al.,  1992)  to  compute  sediment  transport  based  on  grain  size  class. 
These  versions  are  denoted  “HEC-6  Versions”. 

Calculations  for  all  of  the  above  functions  are  made  from  a  common  data  base. 
Sample  input  is  shown  in  Appendix  A.  Thomas  et  al.  have  offered  the  following 
summary  of  the  sediment  transport  equations  and  their  capabilities,  which  is  paraphr^ed 
herein.  While  varying  significantly  in  their  individual  approach,  all  of  the  equations 
share  a  common  functional  form.  This  is  given  as 

GSi  =/(  V,  D,  S,  W,  d,  A ,  Pli ,  s„  T,  Sf) 

where 

V  =  flow  velocity 
D  =  water  depth 
S  =  energy  slope 

W  =  with  of  the  portion  of  the  channel  cross  section  which  is  transporting  bed 
material  sediment 

d,  =  representative  particle  size  for  the  sediment  mixtiue 
dj  =  sediment  particle  size  in  class  I 
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Plj  =  percentage  of  total  sediment  particles  in  size  class  I 
Ss  =  specific  gravity  of  the  sediment  particles 
T  =  temperature  of  the  fluid 
Sf  =  specific  gravity  of  the  fluid 

All  of  the  sediment  transport  fimctions  listed  above  were  originally  supported  in 
the  HEC-6  software  program  (HEC,  1976)  with  the  exception  of  the  Brownlie  equation. 
A  brief  description  of  each  equation  follows. 

Ackers- White  This  is  a  single  grain  size  function  for  sand  bed  streams.  There  is 
also  a  version  which  computes  sediment  transport  by  size  class  using  the  median  grain 
size  for  the  class  as  the  single  grain  size.  The  various  classes  are  then  totaled  to  arrive  at 
an  overall  transport  value.  This  multiple  grain  size  version  is  the  “HEC-6  version”. 

Brownlie  This  is  a  single  grain  size  function  for  sand  transport.  It  is  presented  in 
greater  detail  in  Chapter  1  of  this  report. 

Colby  This  is  a  single  grain  size  function  for  sand  transport  in  streams  and  small 
rivers.  It,  too,  is  covered  in  more  detail  in  Chapter  1. 

Laursen  (Madden)  This  is  a  multiple  grain  size  function  for  sand  bed  transport. 
It  has  also  been  used  for  mixtures  of  sand  and  gravel. 

Laursen  (Copeland)  This  is  a  modification  of  the  Laursen  (Madden)  equation. 
Meyer-Peter  and  Muller  This  is  a  multiple  grain  size  function  for  gravel  bed 

rivers. 

Toffaleti  This  is  a  multiple  grain  size  function  for  sand  bed  rivers. 

Yang  This  is  a  single  size  function  for  sand  transport  in  streams  and  small  rivers. 
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See  Chapter  1  for  a  more  detailed  explanation. 

Quoting  from  the  SAM  user’s  manual  (Thomas  et  al.,  1992),  “This  is  not  an 
exhaustive  list  of  transport  functions.  They  were  selected  from  the  literature  based  on 
experience.  It  does  not  imply  that  those  not  selected  are  deficient.  The  objective  was  to 
provide  designers  with  a’few’  acceptable  methods  which  could  be  supported  for  their  use. 
The  criteria  for  selection  are 

a.  To  cover  a  broad  range  of  particle  sizes. 

b.  To  cover  a  broad  range  of  hydraulic  conditions. 

c.  To  calculate  sediment  transport  by  partitioning  the  mixture  into  size  classes 
and  summing  the  rate  of  each  to  get  the  total  except  when  djo  functions  are  requested. 

d.  To  have  a  history  of  being  reliable  when  used  within  the  range  of  data  for 
which  each  was  calibrated.” 

Together,  this  group  of  functions  represents  a  very  comprehensive 
collection  for  the  prediction  of  sediment  transport,  especially  for  sand  bed  streams.  After 
careful  analysis  of  the  initial  results  from  the  SAM  computer  runs,  10  of  the  13  transport 
functions  were  used  in  the  analyses.  The  two  Meyer-Peter  and  Muller  functions  and  the 
hybrid  Toffaleti-Meyer-Peter  and  Muller  function  were  not  used.  The  Meyer-Peter  and 
Muller  function  is  basically  intended  for  use  on  streams  with  no  appreciable  suspended 
sediment  load  (V anoni,  1975).  Therefore,  this  function  would  not  be  expected  to  perform 
well  on  the  St.  Francis  River  which  is  predominantly  a  sand  bed  stream  with  a  high 
suspended  sediment  load  and  was  removed  from  further  consideration.  The  Toffaleti- 
Meyer-Peter  and  Muller,  which  uses  the  higher  sediment  transport  value  calculated  by  the 
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two  methods,  should  default  to  the  values  of  the  Toffaleti  function  for  most  of  the 
St.  Francis  calculations.  Since  the  Toffaleti  transport  function  is  included  in  the  analysis, 
this  hybrid  function  was  also  not  used.  The  Toffaleti-Schoklitsch  transport  function  also 
combines  a  function  intended  for  high  suspended  loads  with  that  of  one  intended  for  bed 
load  transport  (Vannoni,  1975;  Shultis,  1935).  This  was  the  one  hybrid  function  retained, 
just  for  comparison  to  the  basic  Toffaleti  function. 

An  input  file  for  the  SAM  software  was  prepared  that  represented  all  the  Fisk 
data.  Most  of  the  functions  actually  call  for  the  slope  of  the  energy  grade  line  as  an  input 
parameter.  However,  for  a  river  such  as  the  St.  Francis,  the  assumption  that  the  water 
surface  slope  and  the  slope  of  the  energy  grade  line  are  equivalent  is  valid  (Chow,  1959). 
This  is  so  due  to  the  relatively  gentle  slopes  and  low  energies  involved.  Therefore,  the 
previously  computed  water  surface  slopes  were  input  for  the  slope  values.  It  was  also 
necessary  to  convert  the  temperature  values  from  degrees  Celsius  to  degrees  Fahrenheit. 
This  had  been  done  with  the  SAS  software  and  these  values  were  input.  Bed  material 
data  were  input  as  the  discrete  values  shown  in  Appendix  A.  From  these,  the  SAM 
software  developed  the  grain  size  distribution  curves  and  extracted  the  necessary  grain 
size  information  for  each  transport  function.  An  example  of  the  input  to  and  output  from 
the  SAM  software  is  included  in  Appendix  A. 

The  SAM  software  was  then  used  to  calculate  sediment  transport  values  using  the 
10  selected  sediment  transport  functions.  Values  of  bed  material  sediment  transport  are 
presented  in  tons  per  day.  This  is  consistent  with  the  values  reported  in  the  USGS  field 
measurements.  These  values  are  shown  in  Table  3.  The  procedure  was  repeated 
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for  the  Wilhelmina  Cutoff  and  Clark  Comer  Cutoff  stations.  These  results  are  tabulated 
^  Tables  4  and  5,  respectively.  The  following  abbreviations  are  used  for  the  column 
headings. 

OBS  QS  -  Measured  bed  material  discharge,  computed  by  multiplying  the  total 
measured  load  (which  included  bed  material  load  and  wash  load)  by  1 
minus  the  percent  of  suspended  sediment  smaller  than  0.062  mm 
Material  smaller  than  this  is  in  the  silt-clay  class  of  particle  sizes  which 
largely  constitute  the  wash  load  and  non-cohesive  transport  functions  are 
not  applicable  to  these  size  classes  of  particles. 

TOFF  -  Sediment  transport  calculated  with  the  Toffaleti  transport  function. 
YNG6  -  WES  modified  version  ofthe  Yang  sediment  transport  function.  It 
calculates  sediment  fiansport  by  size  class  and  sums  the  values  for  size 
classes  present. 

AWH6  -  WES  modified  version  of  the  Ackers- White  sediment  transport 
function.  It  calculates  sediment  transport  by  size  class  and  sums  the 
values  for  size  classes  present. 

CBY6  -  WT:S  modified  version  of  the  Colby  sediment  transport  function.  It 
calculates  sediment  transport  by  size  class  and  sums  the  values  for  size 
classes  present. 

TFSH  -  A  hybrid  sediment  transport  function  developed  by  the  WES.  It 
calculates  sediment  transport  by  both  the  Toffaleti  and  Schoklitsch 
methods  and  uses  the  greater  value. 
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BRN  -  Sediment  transport  calculated  with  the  Brownlie  transport  function. 

LSMD  -  Sediment  transport  calculated  with  the  Laursen-Madden  transport 
function. 

LSCP  -  Sediment  transport  calculated  with  the  Laursen-Copeland  transport 
function. 

YG50  -  Sediment  transport  calculated  with  the  Yang  single  grain  size 
transport  function. 

AW50  -  Sediment  transport  calculated  with  the  Ackers- White  single  grain  size 
transport  function. 
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Observed  Water  and  Sediment  Discharge  and  Computed  Bed  Material  sediment 
Discharge  using  10  Different  Transport  Functions  at  the  Fisk,  Mo.  Sampling 
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Table  3  (Continued) 


Water  Flow  OBS  Qs  TOFF 


3500 

45.4 

632 

702 

841 

535 

967 

335 

685 

668 

956 

29 

71 

71 

28 

34 

52 

54 

52 

70 

55 

252 

0.73 

1 

6 

3 

0 

1 

2 

2 

5 

6 

2 

3950 

77.64 

555 

661 

2663 

627 

741 

524 

1531 

495 

644 

3830 

148.91 

640 

813 

1409 

549 

863 

624 

1520 

530 

773 

883 

7980 

167926 

1806 

8834 

2125 

2178 

2417 

1130 

37.4 

54 

29 

132 

21 

57 

0 

635 

254 

0 

0 

4790 

88.72 

447 

387 

746 

285 

548 

0 

3357 

1116 

0 

0 

7410 

270.10 

1432 

1127 

1857 

1665 

3977 

2598 

mm 

2910 

26.71 

364 

424 

1041 

257 

475 

377 

718 

282 

476 

374 

3320 

377.56 

646 

704 

1136 

330 

872 

528 

861 

286 

698 

732 

3630 

104.28 

487 

740 

2464 

■eehbs 

524 

1286 

405 

703 

705 

4410 

599.72 

381 

715 

900 

408 

545 

617 

992 

413 

773 

579 

2510 

104.1 

233 

325 

2119 

270 

314 

328 

1095 

881 

350 

399 

1710 

14.22 

137 

217 

384 

213 

183 

167 

537 

280 

210 

309 

3210 

75.49 

293 

434 

1358 

332 

888 

331 

433 

613 

1960 

21.43 

177 

214 

284 

171 

211 

193 

283 

197 

205 

218 

440 

7.78 

7 

18 

17 

9 

7 

15 

24 

28 

19 

15 

2920 

276.73 

539 

545 

521 

463 

432 

255 

555 

431 

57 


Table  3  ( 

Concluc 

led) 

Water  Flow 

OBSQs 

TOFF 

YNG6 

AWH6 

CBY6 

TFSH 

BRN 

LSMD 

LSCP 

YG50 

AW50 

608 

4.60 

— 

20 

38 

40 

33 

20 

35 

54 

59 

39 

^7 

86 

1 

0 

0 

0 

0 

0 

0 

1 

85 

2.22 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

60 

0.01 

1 

1 

0 

11 

1 

2 

10 

18 

1 

0 

1010 

4.4 

44 

80 

92 

80 

49 

67 

131 

121 

70 

66 

4260 

^  72.35 

518 

1  608 

1261 

581 

658 

563 

1584 

619 

n  609 

641 

2070 

36.89 

569 

314 

17007 

59 

625 

0 

13566 

4436 

0 

0 

2200 

25.19 

299 

314 

1811 

230 

361 

290 

1150 

436 

310 

386 

3150 

42.53 

751 

583 

1348 

449 

888 

537 

1056 

453 

575 

614 

310 

1 

7 

4 

0 

1 

4 

2 

8 

8 

4 

194 

0.53 

0 

3 

1 

0 

HBi 

1 

0 

3 

3 

1 

1 

(Sheei 

t3of  3) 
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Table  4 

Observed  Water  and  Sediment  Discharge  and  Computed  Bed  Material  Sediment 
Discharge  using  10  Different  Transport  Functions  for  the  Wilhelmina  Cutoff,  Mo. 
Sampling  Station 


FLOW 


1470 


TOFF  YNG6 


88.73  1455 


4981 


G6  AWH6 


22002 


33981  6134644 


CBY6  TFSH  BRN  LSMD  LSCP  YG50 


1474  6326  1629 


45651  138361  95021  57373]  10478]  32032]  90340 


1096  887 


180.71 


2030  667.58]  132 


2720  433.00  2263 


6920  1068.00 


434  572 


88  175 


324  323 


1256  2008 


2528 


4375  1299697 


13765  6445910 


2310 

■SI 

12668 

6460 

13658 

9223 

4521 

2290 

351 

122 

5556 

153 

9 

237 

103 

80 

46 

926 

334 

5207 

2370 

1006 

382 

3804 

1565 

7124 

4318 

2592 

1583 

71 

3 

309 

110 

47 

3 

65 

'  25 

31 

0 

254 

94 

2622 

1302 

191 

1  56 

1491 

596 

513  1912 


430  1010 


383 


323  1779 


4317  2209  5961 


2 


112 

47 

477 

242 

136 

65 

315 

209 

865 

2285 

2558 

5085 

33614 

10553 

49374 

19906 

403569 

22654 

75242 

5091 

6662 

314 

206 

7179 

10170 

88 

30 

209 

162 

105 

67 

755 

712 

1369 

2092 

5030 

7450 

7948 

10912 

3288 

3692 

43 

16 

265 

178 

37 

9 

85 

46 

20 

6 

309 

236 

3653 

6217 

244 

199 

1612 

2108 

1292 

1414 

950 

781 

5404 

6596 

602 

378 

722 

684 

4330 

4906 

161 
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Table  4  (Continued) 


FLOW 


202 


379 


6700 


4550 


1240 


316 


1950 


1140 


2370 


338 


1760 


9820 


3460 


1550 


1910 


1680 


668 1  33.001  466 


1159 

416 

886 

813 

Table  4  (Concluded) 

FLOW 

QS 

TOFF 

YNG6 

AWH6 

TFSH 

BRN 

LSMD 

LSCP 

YG50 

AW50 

228 

0.6 

49 

197 

317 

12 

186 

64 

231 

75 

181 

113 

213 

1.20 

23 

67 

43 

25 

33 

35 

46 

74 

69 

43 

800 

60.0 

189 

1327 

22423 

18 

696 

463 

3549 

1560 

3373 

77141 

750 

2718 

1090 

2928 

5259 

7530 

367.0 

3149 

9892 

2621164 

1391 

5626 

3155 

27112 

4774 

7673 

21786 

2250 

806.00 

1655 

4305 

63926 

1125 

3229 

1605 

8825 

2430 

4323 

7141 

2190 

429.0 

1226 

3236 

16994 

621 

2643 

1304 

5090 

1370 

3258 

3632 

6051 

19063 

1935 

5072 

3054 

5975 

7605 

603 

26.00 

166 

253 

230 

106 

260 

239 

199 

253 

212 

289 

1.77 

229 

132 

36 

239 

72 

79 

85 

230 

112 

(Sheet  3  of  3) 
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Table  5 

Observed  Water  and  Sediment  Discharge  and  Computed  Bed  Material  Sediment 
Discharge  using  10  Different  Transport  Functions  for  the  Clark  Corner  Cutoff,  AR. 
Sampling  Station 

FLOW  I  QS  TOFF  |yNG6  AWH6  CBY6  TFSH~  f  RPM  1 1  CMn  ll  [77 

200  0 


BRN 

LSMD 

1 

0 

78 

1 

196 

36 

YG50 

AW50 

3 

1 

37 

49 

243 

145 

25 

12 

640 

1705 

1 

1  0 

6266  38460 


84431  67561 


14 

3 

152 

80 

58 

69 

82 

89 

91 

100 

21 

19 

3 

2 

8 

3 

19 

10 

20 

14 

36 

46 

3216 

6111 

369 

507 

109 

58 

12 

5 

46 

13 

29 

6 

77 

62 

1 

0 

■K^l 

mm\ 

3521 1 

4501 1 

(Continued 
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Table  5  (Concluded) 


FLOW 

QS 

2190 

26 

26300 

3544 

mmm 

1110 

32 

13300 

2986 

10300 

2102 

5610 

42 

10000 

21300 

749 

2290 

12 

9230 

2213 

21400 

3494 

15900i 

440 

15800 

2399 

12200 

1069 

3950 

362 

3260 

269 

7920 

887 

8140 

697 

9190 

1533 

6470 

740 

16200 

■BEII 

1705 

45 

3600 

24 

718 

5 

2540 

57 

8770 

2253 

3230 

836 

3 

1180 

8 

2990 

105 

2410 

37 

15300 

3397 

10100 

891 

7900 

667 

6840 

279 

2660 

67 

397 

7.8 

YG50  lAWSO 


2898 

210000 

2028 

1  100000 

523 

541 

59 

89 

1 

1  1 

The  output  in  Tables  3-5  were  used  as  a  basis  for  further  research  as  discussed  in 
Chapter  4. 
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4  Analysis  and  Results 


Sediment  Transport 


In  order  to  evaluate  the  field  data,  a  method  had  to  be  devised  to  measure  the 
value  of  the  data,  both  as  a  complete  set  and  incrementally.  The  sediment  transport  values 
calculated  with  the  SAM  software  and  depicted  in  Tables  3-5  were  generated  to  provide  a 
means  for  assisting  in  answering  the  question,  “How  much  sediment  data  are  necessary?” 
While  it  is  customaiy  to  use  the  field  data  to  select  a  sediment  transport  function  (Mien, 
1995).  It  is  also  possible  that  the  transport  functions  could  be  used  to  evaluate  the 
incremental  and  total  value  of  the  data  sets. 

In  order  to  accomplish  this,  it  was  first  necessary  to  determine  which  transport 
functions  best  reproduced  the  measured  field  data.  Water  discharge  values  were  plotted 
versus  values  from  the  10  selected  transport  functions  and  the  observed  data.  For  ease  of 
reading,  these  were  plotted  five  to  a  graph.  Figure  6  depicts  a  plot  of  water  discharge  vs. 
sediment  discharge  for  the  observed  sediment  discharge  and  five  predictive  equations  for 
the  Fisk,  MO,  Station.  This  shows  graphically  how  well  the  equations  reproduced  the 
observed  data.  Figure  7  depicts  similar  relationships  for  the  other  five  predictive 
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Figure  6.  Water  vs.  sediment  discharge,  five  transport  functions,  Fisk,  MO 
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Figure  7.  Water  vs.  sediment  discharge,  other  transport  functions,  Fisk,  MO 


equations.  Similarly,  Figures  8  and  9  depict  the  water  discharge  vs.  sediment  discharge 
relations  for  five  predictors  for  the  Wilhelmina  Cutoff  station.  Figures  10  and  1 1  show 
the  predictors  for  the  Clark  Comer  station.  These  visual  plots,  while  informative,  do  not 
present  results  in  a  quantitative  manner  which  establishes  a  firm  basis  for  evaluation  of 
which  of  the  functions  best  fit  the  observed  sediment  data.  Other  quantitative  methods 
were  applied  to  this  end  and  the  results  are  discussed  below. 


Regression  Analysis 


While  visual  observation  of  the  water  vs.  sediment  discharge  plots  showed  general 
comparisons  of  the  measured  vs.  computed  sediment  transport,  it  was  insufficient  to  rank 
the  performance  of  the  various  sediment  transport  prediction  functions.  To  this  end,  a 
simple  regression  analysis  was  made  for  each  transport  function  as  compared  to  the 
measured  bed  material  load.  The  correlation  coefficients,  R,  were  computed  for  each  set 
of  measured  vs.  computed  data.  Using  these,  the  coefficient  of  determination,  R^,  were 
computed  in  each  case.  The  coefficient  of  determination  is  the  value  of  the  correlation 
coefficient,  squared.  This  value,  R^,  is  a  ratio  of  the  explained  variation  in  the  dependent 
variable,  in  this  case  the  computed  sediment  discharge,  to  the  total  variation  of  the 
dependent  variable.  The  transport  functions  were  then  ranked  firom  one  to  ten,  with  “one” 
indicating  the  highest  coefficient  of  determination  and  “ten”  the  worst.  These  results  are 
shown  in  Table  6  for  Fisk,  MO. 
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Figure  8.  Water  vs.  sediment  discharge,  five  transport  functions,  Wilhelmina  MO 
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Figure  9.  Water  vs.  sediment  discharge,  other  transport  functions,  Wilhelmina,  MO 
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Table  6 

Coefficients  of  Determination  for  Each  Transport  Function  Regressed  Against 
the  Measured  Values  Fisk,  MO 


Measured  Sediment  Discharge 
Regressed  Against 

Coefficient  of  Determination,  R^ 

Rank 

Toffaleti 

32.16 

5 

Yang-HEC6 

49.92 

1 

Ackers-White-HEC6 

35.56 

4 

Colby-HEC6 

25.27 

7 

Toffaleti-Schoklitsch 

39,75 

2 

Brownlie 

36.00 

3 

Laursen-Madden 

27.49 

6 

Laursen-Copeland 

15.44 

8 

Yang,  ds,, 

6.86 

9 

Ackers-White,  djo 

0.09 

10 

It  can  be  seen  that  even  the  top  rated  transport  function,  Yang-HEC6,  has  less  than 
one  half  of  its  variation  explained.  This  indicates  that  it  is  very  poorly  correlated  with  the 
observed  data.  The  Yang-HEC6  function  was,  however,  clearly  better  than  any  of  the 
others.  That  it  performed  best  is  not  surprising  since  it  was  developed  for  sand-bed 
streams.  The  WES  modifications  to  allow  sediment  transport  calculations  by  size  class 
greatly  improve  the  function  in  this  case.  The  unmodified  single  grain  size  version 
finished  rated  ninth.  The  Ackers- White  single  grain  size  function  was  the  worst 
performer,  with  almost  no  explained  variation.  It  therefore  has  almost  no  correlation  to 
the  observed  data.  This  equation  has  previously  been  reported  to  be  a  poor  predictor  for 
broadly  graded  samples  (Prasuhn  et  al.,  1987).  In  addition,  it  was  developed  for  a 
maximum  grain  size  of  2.5  mm,  which  was  exceeded  by  a  number  of  the  samples  at  Fisk. 
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The  Yang  single  grain  size  function  was  also  a  very  poor  performer.  Careful  study  of  the 
data  (see  Appendix  A)  revealed  that  a  wide  variability  existed  in  the  bed  material 
samples.  The  maximum  size  ranged  from  16  mm  to  .0625  mm.  The  average  djo  size  was 
0.3549  mm.  However,  there  were  five  samples  where  the  maximum  size  present  in  the 
bed  was  16  mm,  three  samples  where  the  maximum  size  was  8mm  and  three  samples 
where  the  maximum  size  was  4  mm.  This  wide  variability  of  bed  material  size  limits  the 
ability  of  any  one  sediment  transport  function  to  be  a  good  predictor,  especially  the  single 
grain  size  functions.  In  fact  it  is  quite  remarkable  that  the  Brownlie  single  grain  size 
function  did  as  well  as  it  did  (third  overall).  The  average  djo  for  all  the  observations  is 
classified  as  a  fine  sand.  However,  the  maximum  djo  is  6.0  mm.  This  is  in  the  range  of 
fine  gravel.  Therefore,  the  wide  variability  in  bed  material  size  limits  the  ability  of  any  of 
the  transport  functions  to  accurately  reproduce  the  observed  data.  With  this  limitation  in 
mind,  the  results  were  kept  and  the  analysis  continued.  If  a  procedure  could  be  developed 
which  covered  this  eventuality  as  well  as  more  conventional  data,  then  it  would  be  robust 
indeed.  More  conventional  results  were  obtained  for  Wilhelmina  and  Clark  Comer 
Cutoffs,  as  shown  in  Tables  7  and  8. 

These  sediment  data  were  far  less  variable  than  those  at  Fisk.  The  average  djo  was 
0.1775  mm.  The  maximum  grain  size  was  2.0  mm  for  two  samples.  The  remainder  of  the 
samples  had  a  maximum  grain  size  of  1 .00  mm.  This  lack  of  variability  is  reflected  in  the 
much  higher  coefficient  of  determination  values  for  Wilhelmina  as  compared  to  those  for 
Fisk.  It  should  be  noted  that  the  Ackers- White  function  continued  to  suffer  from  its 
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Table  7 

Coefficients  of  Determination  for  £ach  Transport  Function  Regressed  Against 
the  Measured  Values  Wilhelmina  Cutoff,  MO 


Measured  Sediment  Discharge 
Regressed  Against 

Coefficient  of  Determination, 

Rank 

Toffaletl 

62.97 

5 

Yang-HEC6 

65.73 

4 

Ackers-White-HEC6 

8.24 

10 

Colby-HEC6 

68.53 

2 

Toffaleti-Schoklitsch 

62.91 

6 

Brownlie 

72.82 

1 

Laursen-Madden 

61.62 

7 

Laursen*CopeIand 

60.76 

8' 

Yang,  djo 

65.75 

3 

Ackers-White,  dso 

10.85 

9 

Table  8 

Coefficients  of  Determination  for  Each  Transport  Function  Regressed  Against 
the  Measured  Values  Clark  Comer  Cutoff,  AR 

{Measured  Sediment  Discharge 
Regressed  Against 

Coefficient  of  Determination,  R* 

Rank 

Toffaleti 

78.44 

2 

Yang-HEC6 

68.97 

4 

Ackers-Whlte-HEC6 

22.79 

9 

Colby-HEC6 

52.33 

8 

Toffaleti-Schoklitsch 

78.46 

1 

Brownlie 

71.69 

3 

Laursen-Madden 

65.87 

6 

Laursen-Copeland 

68.57 

5 

Yang, 

61.98 

7 

Ackers-White,  dso 

5.47 

10 
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sensitivity  to  broadly  graded  bed  material  in  both  its  single  and  multiple  grain  size 
versions.  The  Brownlie  function  finished  first  at  Wilhelmina.  This,  coupled  with  its 
third  place  finish  at  Fisk  indicate  that  it  is  a  versatile  function  and  indicate  the  strength  of 
the  regression-based  approach  used  to  develop  this  function  (Brownlie,  1981). 

The  highest  coefficients  of  determination  of  the  three  stations  analyzed  were 
computed  for  Clark  Comer  Cutoff.  The  variation  in  bed  material  was  limited  at  this 
station,  and  this  undoubtedly  contributed  to  the  overall  better  performance  of  the 
transport  functions.  In  these  data,  one  sample  had  a  maximum  grain  size  of  8.0  mm  and 
seven  samples  had  a  maximum  size  of  2.0  mm.  The  rest  were  fairly  uniform.  Even  the 
2.0  mm  value  is  in  the  medium  to  coarse  sand  range.  The  average  djofor  this  station  was 
0.2748  mm.  Therefore,  the  selected  transport  functions  should,  and  did,  by  and  large 
perform  quite  well  at  this  station.  The  hybrid  Toffaleti-Schoklitsch  function  was  the  top 
rated  function.  It  should  be  noted  that  it  is  only  very  marginally  better  than  the  Toffaleti 
transport  function.  This  indicated  that  the  Toffaleti  function  was  used  for  virtually  all  of 
the  computations,  as  would  be  expected.  Brownlie  also  did  quite  well  again  and  finished 
third.  The  two  Ackers- White  functions,  as  was  the  case  in  general  for  all  the  stations, 
were  very  poor  performers.  This  is  all  the  more  puzzling  since,  based  on  the  conditions 
for  which  it  is  recommended,  it  should  have  been  a  good  predictor.  The  likely 
explanation  lies  in  the  fact  that  the  Ackers- White  function  is  based  on  a  mobility  number 
and  a  dimensionless  grain  diameter.  The  mobility  number  is  a  stream  power  derivative 
and  should  function  as  well  as  the  Yang  stream  power  based  function.  The  dimensionless 
grain  diameter,  however,  seems  more  problematic.  In  comparing  the  function  to  actual 
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field  data.  Ackers  and  White  had  to  abandon  the  d5o  grain  size  in  favor  of  the  d35  grain 
size  in  order  to  best  match  the  selected  field  data.  The  SAM  program  uses  the  d35  grain 
size.  However,  for  the  St.  Francis  data,  this  value  may  be  quite  variable.  As  this  was  the 
only  function  that  used  the  d35  value  and  it  was  a  poor  predictor  in  general,  this  possibility 
was  not  investigated  further. 

The  above  analyses  were  performed  using  the  full  eleven  years  of  record  at  each 
station.  It  was  recognized  that  to  meet  the  research  objective,  the  minimum  number  of 
years  of  record  required  to  select  the  best  transport  function  must  be  determined.  To  this 
end,  the  regression  analyses  were  repeated  for  each  station.  However,  the  analysis  at  each 
station  was  broken  into  eleven  segments.  The  first  segment  contained  only  the  first  year 
of  record.  Segment  number  two  contained  the  first  two  years  of  record  and  so  on  until 
segment  number  eleven  contained  the  full  record;  This  was  repeated  for  each  of  the  three 
stations.  The  sediment  transport  equations  were  ranked  for  each  segment  as  described 
above.  The  notion  here  is  that  if,  for  instance,  the  beginning  order  of  rank  of  the  transport 
functions  does  not  change  as  additional  data  are  added,  then  the  additional  data  did  not 
aid  in  making  a  better  selection  of  a  transport  function  that  represents  the  data.  The  other 
extreme  possibility  was  that  the  order  would  continuously  shift,  indicating  no  pattern.  It 
was  also  necessary  to  select  the  best  transport  function  for  use  in  application  of  the 
entropy  principles  described  below. 

The  coefficients  of  determination  for  the  Fisk,  MO,  station  for  each  segment  of 
data  are  tabulated  in  Table  9.  (Note.  Abbreviations  used  in  Tables  9-14  are  the  same  as 
those  previously  defined  for  Tables  3-5.)  They  are  plotted  in  Figure  12  (best  five 
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predictors)  and  Figure  13  (other  five  predictors)  versus  the  number  of  years  of  data. 

Table  10  shows  the  rank  order  of  each  transport  function  as  they  varied  with  an  increasing 
number  of  data  points  in  the  data  set.  The  top  four  functions  are  high-lighted  in  bold  for 
each  segment.  These  data  are  plotted  in  Figure  14  (best  five  predictors)  and  Figure  15 
(other  five  predictors). 


Table  9 

Coefficient  of  Determination  Between  Measured  and  Observed  Sediment 
Transport  with  Increasing  Number  of  Years  of  Data,  Fisk,  MO 


YNG6 

AWH6 

CBY6 

LSMD 

AW50 

1 

48.64 

51.28 

45.70 

63.79 

50.53 

55.90 

49.38 

52.61 

62.36 

45.42 

2 

45.08 

57.23 

36.88 

34.22 

44.18 

45.53 

52.15 

44.23 

2.22 

3 

47.08 

58.22 

37.31 

47.00 

46.10 

52.53 

46.64 

3.89 

4 

57.80 

52.67 

44.22 

45,98 

52.31 

3.16 

0.24 

5 

36.76 

42.09 

3.88 

0.16 

6 

35.67 

25.98 

35.54 

43.93 

4.28 

0,15 

7 

34.18 

49.40 

36.85 

25.00 

41.26 

35.37 

41.42 

0.13 

8 

38.62 

25.99 

41.31 

35.16 

‘  0.15 

9 

33.36 

49.38 

35.37 

25.91 

40.30 

35.65 

40.03 

31.41 

6.09 

0,16 

10 

34.26 

50.14 

35.58 

26.69 

41.21 

36.58 

32.12 

6.58 

0.12 

11 

32.16 

49.92 

35.56 

36,00 

■OS 

6.86 

0.09 

Table  10 

Relative  Rank  of  Sediment  Transport  Equations  with  Increasing  Number  of 
Years  of  Data,  Fisk,  MO 
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Coefficient  of  Determination  Ualue 


Number  of  Years  in  Data  Set 


Figure  14.  Number  of  years  of  data  vs.  rank  of  the  best  five  transport  functions 
Fisk,  MO 


Rank  of  Function 


It  can  be  seen  that  selection  of  a  transport  function  based  on  the  top  ranked  function  after 
one  year  of  data  (Colby-HEC6)  would  result  in  a  very  poor  long-term  selection.  This  is 
proven  as  it  drops  to  ninth  in  year  two  and  never  climbs  higher  than  seven  during  the 
period  of  record.  Selection  of  the  top  rated  function  using  two  years  of  data  would  result 
in  the  same  choice  as  that  based  on  eleven  years  of  data.  Note  that  after  seven  years, 
there  is  also  less  shifting  of  rank  among  the  top  four  ranked  functions. 

The  coefficients  of  determination  for  the  Wilhelmina  Cutoff  Station  for  each 
segment  of  data  are  tabulated  in  Table  1 1 .  They  are  plotted  in  Figure  16  (best  five 
predictors)  and  Figure  17  (other  five  predictors).  Table  12  shows  the  rank  order  of  each 
transport  function  as  they  varied  with  an  increasing  number  of  data  points  in  the  data  set, 
again  with  the  top  four  highlighted  for  each  segment  These  data  are  plotted  in  Figure  18 
(best  five  predictors)  and  Figure  19  (other  five  predictors). 

The  choice  of  the  top  rated  function  based  on  one  year  of  data,  as  was  the  case  at 
Fisk,  would  not  result  in  the  best  long-term  selection.  Selection  of  the  top  rated  function 
using  two  years  of  data  would  result  in  the  same  choice  as  that  based  on  eleven  years  of 
data,  which  was  the  same  conclusion  drawn  from  the  Fisk  data.  Also  as  at  Fisk,  note  that 
after  seven  years  there  is  much  less  shifting  of  rank  among  the  top  four  ranked  functions. 

The  coefficients  of  determination  for  the  Clark  Comer  Cutoff  Station  for  each 
segment  of  data  are  tabulated  in  Table  13.  They  are  plotted  in  Figure  20  (best  five 
predictors)  and  Figure  21  (other  five  predictors).  Table  14  shows  the  rank  order  of  each 
transport  function  as  they  varied  with  an  increasing  number  of  data  points  in  the  data  set. 
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Coefficient  of  Determination  Between  Measured  and  Observed  Sediment 


YEAR 

TOFF 

YGH6 

AWH6 

CBY6 

TFSH 

BRN 

LSMD 

LSCP 

YG50 

AW50 

1 

21.56 

14.64 

7.24 

10.27 

20.12 

15.46 

8.85 

2 

34.09 

31.43 

4.29 

27.06 

38.66 

20.90 

27.80 

1.42 

3 

35.33 

25.35 

40.78 

44.75 

29.97 

2.89 

4 

35.65 

25.33 

40.70 

32.31 

27.05 

33.48 

3.93 

5 

9.55 

24.82 

40.10 

44.17 

33.31 

28.45 

34.21 

4.45 

6 

45.85 

45.99 

8.08 

40.73 

47.95 

51.85 

42.28 

36.87 

43.93 

5.96 

7 

44.41 

45.74 

8.22 

39.32 

46.77 

51.17 

41.82 

36.56 

43.63 

6.11 

8 

63.88 

66.01 

8.26 

68.52 

63.89 

73.13 

62.33 

61.38 

65.93 

15.34 

9 

63.54 

65.86 

8.44 

63.35 

63.54 

72.89 

62.34 

61.31 

65.74 

15.51 

10 

63.79 

62.26 

8.80 

68.67 

63.72 

73.20 

62.81 

61.65 

66.13 

15.97 

11 

62.97 

65.73  8.24 

68.53 

62.91 

10.85 

Table  12 

Relative  Rank  of  Sediment  Transport  Equations  with  Increasing  Number 
of  Years  of  Data,  Wilhelmina  Cutoff,  MO 
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Table  13 


Coefficient  of  Determination  Between  Measured  and  Observed  Sediment  Transport 

_ X _ •  'fc.T  «  MW.  _  _  _  _  * 


TOFF  YNG6  AWH6  CBY6  TFSH  BRN  LSMD 

0.62 

0.23 1 

0.74| 

11.32 

0.63 1 

0.04 1 

0.79 

95.58 

94.31 

88.78 

93.831 

95.821 

94.98 

89.02 

88.25 

94.07 

95.13 

89.33 

88.59 

94.21 

88.80 

86.90 

95.81 

60.76 

88.94 

82.00 

93.91 

87.23 

80.61 

92.41^ 

52.79 

87.43 

76.92 

91.95 

86.68 

77.73 

92.09 

52.56 

86.82 

76.33 

91.60 

86.16 

74.88 

25.23 

52.60 1 

86.30 

73.00 

84.71 

85.90 

86.07 

73.19 

84.69 

85.61 

74.82 

25.21 

52.31 

85.78 

73.08 

84.43 

78.44 

68.97 

22.79 

52.33 

78.46 

71.69 

65.87 

ZE9j 


AW50  I 


0.41 


83.89 


84.41 


84.79 


84.75 


83.15 


Table  14 

Relative  Rank  of  Sediment  Transport  Equations  with  Increasing  Number  of 
Years  of  Data,  Clark  Corner  Cutoff,  AR  _ 

YEAR  ItOFF  |yNG6  |aWH6  |cBY6  jxFSH  [bRN  IlSMD  IlSCP  IyGSO  l^( 


YG50  AW50 
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Figiire  16.  Variability  of  the  coefficient  of  determination  values  as  more  data  are  added, 
best  five  transport  functions,  Wilhelmina,  MO 


Coefficient  of  Determination  Value 


Rank  of  Function 


Rank  of  Function 


Figure  19.  Number  of  years  of  data  vs.  rank  of  the  other  transport  functions, 
Wilhelmina,  MO 


The  top  four  are  again  highlighted  in  bold  for  each  segment.  These  data  are  plotted  in 
Figure  22  (best  five  predictors)  and  Figure  23  (other  five  predictors). 

As  ■was  true  at  Fisk  and  Wilhelmina,  selection  of  the  top  rated  function  based  on 
one  year  of  data  would  result  in  a  poor  long  term  selection.  Selection  of  the  top  rated 
function  using  two  years  of  data  would  also  result  in  a  poor  long-term  selection.  This 
result  is  quite  different  from  the  conclusion  drawn  from  the  previous  two  stations.  That 
is,  that  the  selection  based  on  two  years  of  data  results  in  the  same  choice  as  that  based  on 
eleven  years  of  data.  In  fact,  it  is  only  after  eight  years  that  the  selection  of  the  top  rated 
function  results  in  the  best  long  term  selection.  It  is,  however,  true  that  the  variability  of 
the  top  four  rated  functions  is  minimal  after  eight  years,  same  as  for  the  previous  two 
stations. 

For  this  analysis  it  is  assumed  that  the  top  ranked  transport  function  based  on  the 
full  eleven  years  of  record  is  the  long  term  best  choice  with  which  to  make  predictions  of 
sediment  transport  at  a  given  station.  This  is  logical  since  such  a  function  has  endured 
the  rigors  of  variability  inherent  in  sedimentation  data.  Indeed,  the  accepted  and  usual 
method  of  selecting  a  transport  function  is  to  compare  several  of  them  to  observed  data 
by  some  method  and  select  the  one  that  best  reproduces  the  observed  data  (Vanoni,  1975, 
USAGE,  1989,  Julien,  1995).  However,  as  can  be  seen  from  the  preceding  analysis,  a 
different  selection  would  be  made  at  the  same  station  depending  on  how  much  data  were 
available  for  inclusion  in  the  selection  process.  At  Fisk  and  Wilhelmina  two  years  of  data 
are  sufficient  to  select  the  best  transport  function  based  on  the  full  eleven  years  of  record. 
Eight  years  of  data  are  necessary  at  Clark  Comer  to  select  the  best  function.  Variability 


91 


Figure  22.  Number  of  years  of  data  vs.  rank  of  the  best  five  transport  functions,  Clark 
Comer,  AR 
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Rank  of  Function 


among  the  top  four  ranked  functions  was  less  after  seven  to  eight  years. 

So  then,  referring  back  to  the  research  objective,  i.e.,  to  determine  how  much  data 
are  necessary ,  the  answer  is  not  entirely  clear.  A  range  can,  however,  be  identified  which 
does  answer,  in  broad  terms,  the  objective  question.  This  range  is  fi-om  two  to  eight  years 
for  the  St.  Francis  river  data  set. 

This  would  indicate  that  if  one  were  to  collect  eight  years  of  sediment  data,  on  a 
monthly  basis,  that  this  would  constitute  a  sufficient  data  set  to  select  the  best  sediment 
transport  equation  for  predicting  sediment  transport  in  a  given  reach  of  the  St.  Francis 
River.  This  statistical  evidence,  however  convincing,  seems  insufficient  to  stand  alone  as 
positive  proof  that  eight  years  of  record  are  all  that  is  necessary. 


Entropy  Calculations 


The  work  of  the  previous  section  defined  a  range  (2-8  years)  which  encompassed 
the  minimum  amount  of  data  necessary  to  determine  the  sediment  transport  function  that 
best  reproduced  the  observed  data.  However,  while  providing  definitive  guidance  that 
has  not  previously  been  available,  this  range  still  seemed  too  broad.  The  value  added 
associated  with  collection  of  the  fully  necessary  data  set  was  also  not  quantified.  In  order 
to  narrow  the  range,  earlier  works  in  the  application  of  the  principle  of  entropy  to  water 
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resources  (see  Chapter  1)  were  used  as  a  basis  to  formulate  an  independent  evaluation  of 
the  amount  of  data  needed  to  provide  the  most  information  for  the  selected  three  stations 
in  the  St.  Francis  Basin.  In  addition,  the  principle  of  maximum  entropy  was  applied  to 
quantify  the  adequacy  of  the  data  sets. 

To  express  the  uncertainty  in  the  computed  sediment  transport  mathematically,  it 
was  assumed  that  the  sediment  transport  is  defined  as  a  random  variable  X  which  has 
values  of  x,,  X2,....x„,  with  probabilities  Pi,P2v-Pn.  P(X=Xi)=p,, 
P(X=X2)=P2v-P(X=Xn)=Pn.  P(x)  is  the  probability  distribution  of  X  satisfying 


n 


(16) 


/=i 


After  Shannon  (1948a,  1948b),  the  entropy  is  defined  as: 


H{_P)=HiJC)  =  -Y,p\np. 

i=\ 


(17) 


where  H(P)  is  the  entropy  of  the  probability  distribution  P=  (Pi,P2,—Pn)-  Using  equation 
(17),  the  incremental  and  total  entropy  could  be  computed.  This  would  reveal  the  rate  of 
information  accumulation  and  the  point  of  maximum  information. 
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Singh  (1992)  has  stated  that  entropy  is  a  measure  of  the  expected  information 
content  of  space-time  measurements.  Since  Jaynes  (1957)  has  also  stated  that  entropy  and 
uncertainty  are  synonymous,  this  connection  is  not  immediately  obvious.  However,  it  has 
been  proven  by  many  authors  (Shannon,  1948;  Jaynes,  1957;  Singh,  Fiorentino,  1992) 
that  when  entropy  is  maximized,  the  resulting  data  set  is  the  least  biased.  Put  another 
way,  as  long  as  entropy  is  increasing,  additional  samples  reduce  the  bias  of  the  population 
and  are  ever  more  valuable  in  that  they  provide  more  information.  Once  the  entropy 
value  stabilizes,  additional  information  is  less  valuable.  This  may  seem  counter  intuitive 
until  one  puts  this  statement  in  perspective.  For  example,  let  us  assume  that  there  are  two 
data  stations,  station  X  and  station  Y.  At  station  X  there  is  perfect  understanding  of  the 
processes  in  which  we  are  interested  in  modeling.  At  station  Y  we  know  nothing  of  the 
processes  in  which  we  are  interested  in  modeling.  Furthermore,  let  us  say  that  we  are  to 
be  allowed  to  take  a  fixed  number  of  additional  data  samples,  Z.  The  question  is,  where 
to  take  them?  Since  we  possess  perfect  knowledge  at  station  X,  there  is  no  value  in 
taking  additional  samples  at  station  X.  However,  data  samples  taken  at  station  Y  would 
be  extremely  valuable,  since  we  have  no  knowledge  here.  Each  additional  data  sample  at 
station  Y  would  increase  our  imderstanding  of  the  processes  there,  up  to  a  point.  That 
point  is  the  point  of  maximum  total  entropy. 

Therefore,  the  data  sets  for  Fisk,  Wilhelmina  and  Clark  Comer  were  analyzed  by 
looking  at  the  incremental  total  entropy  after  each  year  of  data  were  collected.  The  point 
at  which  the  entropy  definitely  ceased  to  increase  would  define  the  minimum  amovmt  of 
sediment  data  necessary  that  provided  the  maximum  amount  of  information  and  least 
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biased  population.  If  the  entropy  continually  increased,  then  the  data  set(s)  would  not 
contain  the  maximum  amount  of  information. 

The  procedure  for  calculating  the  entropy  at  each  station  consisted  of  several 
steps.  First,  the  data  to  be  used  in  the  calculation  were  selected.  These  data  consisted  of 
the  best  calculated  values  for  each  station,  best  being  defined  as  the  transport  function 
that  was  the  best  predictor  based  on  the  regression  analyses  performed  above.  The 
selected  transport  functions  were  Yang-HEC6  for  the  Fisk,  Station,  Brownlie  for  the 
Wilhelmina  station  and  Toffaleti-Schoklitsch  for  the  Clark  Comer  station.  Next,  the  data 
were  presented  as  a  discrete  distribution.  This  was  generated  by  dividing  the  data  into 
intervals.  The  number  of  intervals  selected  was  based  on  trying  to  achieve  a  reasonable 
number  of  intervals.  If  too  few  were  selected,  the  data  would  all  be  clumped  in  these 
intervals  and  the  distribution  might  well  be  biased.  Too  many  intervals  would  lead  to 
uimecessary  calculations  without  changing  the  actual  results.  The  Fisk  station  seemed  to 
offer  no  real  problems  along  these  lines.  The  data  values  fit  within  a  range  of  0-4000. 
Therefore,  forty  intervals  resulted  in  an  interval  width  of  100  (Figures  20-34).  The  other 
two  stations  offered  a  bit  of  problem  in  that  the  range  of  values  were  much  greater.  At 
Wilhelmina,  the  data  values  fit  within  a  range  of  0-26,000  and  at  Clark  Comer  ranged 
from  0-38,000.  Therefore,  the  data  at  these  stations  were  distributed  using  thirty  and 
fifty-two  intervals  at  Wilhelmina  and  forty-five  and  seventy-six  at  Clark  Comer.  The 
different  number  of  intervals  were  used  to  investigate  sensitivity  to  this  decision.  The 
number  of  intervals  changed  the  absolute  value  of  the  entropy,  but  not  the  trends.  The 
higher  number  of  intervals  were  selected  for  analysis  and  are  presented  in  the  figures 
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showing  the  transport  function’s  distribution  at  Wilhelmina  and  Clark  Comer 
(Figures  35-56). 

The  segmented  data  were  plotted  versus  various  probability  distributions.  Several 
distributions  initially  seemed  to  fit  the  data  well,  notably  the  Lognormal,  Gamma  and 
Weibull.  On  closer  inspection,  the  Weibull  distribution  was  selected  as  the  most 
representative  of  the  three  distributions.  Figures  24-34  depict  the  Yang-HEC6  computed 
values  (best  predictor  based  on  regression  analysis)  divided  into  forty  intervals  and  fitted 
to  the  Weibull  distribution  for  the  Fisk,  MO,  Station.  Figures  35-45  depict  the  Brownlie 
computed  values  (best  predictor  based  on  regression  analysis)  divided  into  fifty-two 
intervals  and  fitted  to  the  Weibull  distribution  for  the  Wilhelmina  Cutoff  station. 

Figures  46-56  depict  the  Toffaleti-Schoklitsch  computed  values  (best  predictor  based  on 
regression  analysis)  divided  into  seventy-six  intervals  and  fitted  to  the  Weibull 
distribution  for  the  Clark  Comer  Cutoff  station.  Next,  the  number  of  occurrences  in  each 
interval,  for  each  data  segment  (i.e.,  data  including  one  year,  two  years,  etc.)  were 
tabulated.  These  were  input  into  a  Quattro-Pro  spreadsheet  and  the  probability  of 
occurrence,  p; ,  for  each  interval  calculated  by  dividing  the  number  of  occurrences  by  the 
total  number  of  data  points.  Using  these  probabilities  and  equation  (17),  the  total  entropy 
was  calculated  for  each  data  segment.  In  the  case  of  Wilhelmina  and  Clark  Comer,  two 
sets  of  computations  were  made  using  the  two  different  intervals.  These  entropy  values 
(maximums  boldface)  and  their  variation  with  the  amoimt  of  data  collected  are  shown  in 
Tables  15  and  16. 
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Figure  24.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for 
one  year  of  data,  Fisk,  MO 
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Figxrre  25.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for 
two  years  of  data,  Fisk,  MO 
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Figiire  26.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for 
three  years  of  data,  Fisk,  MO 
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Figure  28.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for 
five  years  of  data,  Fisk,  MO 
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Figure  29.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for 
years  of  data,  Fisk,  MO 
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Figure  3 1 .  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for 
eight  years  of  data,  Fisk,  MO 
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Figure  33.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution  for  ten 
years  of  data,  Fisk,  MO 
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Figure  34.  Yang-HEC6  transport  function  results  fitted  to  the  Weibull  distribution 
eleven  years  of  data,  Fisk,  MO 
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Figure  35.  Brownlie  transport  function  results  fitted  to  the  Weibuli  distribution  for 
year  of  data,  Wilhelmina,  MO 
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Figiire  36.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for  two 
years  of  data,  Wiihelmina,  MO 
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Figure  37.  Brownlie 
years  of  data,  Wilhelr 
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Figure  38.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for  four 
years  of  data,  Wilhelmina,  MO 
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Figure  40.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for  six 
years  of  data,  Wilhelmina,  MO 
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Figure  42.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for  eight 
years  of  data,  Wilhelmina,  MO 
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Figure  43.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for  nine 
years  of  data,  Wilhelmina,  MO 
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Figxire  44.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for  ten 
years  of  data,  Wilhelmina,  MO 
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Figure  45.  Brownlie  transport  function  results  fitted  to  the  Weibull  distribution  for 
eleven  years  of  data,  Wilhelmina,  MO 


120 


TFSH  One  Year  of  Data 


(X  10000) 


Figure  46.  TofFaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  one  year  of  data,  Clark  Comer,  AR 
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Figure  47.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  two  years  of  data,  Clark  Comer,  AR 
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Figure  48.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  three  years  of  data,  Clark  Comer,  AR 
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Figure  49.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  four  years  of  data,  Clark  Comer,  AR 
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Figure  50.  TofFaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  five  years  of  data,  Clark  Comer,  AR 
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Figure  5 1 .  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  six  years  of  data,  Clark  Comer,  AR 
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Figure  52.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  seven  years  of  data,  Clark  Comer,  AR 
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Figure  53.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  eight  years  of  data,  Clark  Comer,  AR 
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Figure  54.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  nine  years  of  data,  Clark  Comer,  AR 
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Figure  55.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  ten  years  of  data,  Clark  Comer,  AR 
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Figure  56.  Toffaleti-Schoklitsch  transport  function  results  fitted  to  the  Weibull 
distribution  for  eleven  years  of  data,  Clark  Comer,  AR 


Table  15 

Calculated  Entropy  for  Each  Station,  by  Years  of  Data  (40  Intervals  for  Fisk, 

30  Intervals  for  Wilhelmina  and  45  Intervals  for  Clark  Comer) 

Number  of  Years  of  Data 

FISK,  MO 

WILHELMINA 

CLARK  CORNER 

1 

1 

0.5 

2 

1.04 

1.27 

0.83 

3 

1.11 

1.3 

0.72 

4 

1.17 

1.28 

0.56 

5 

1.32 

1.28 

0.56 

6 

1.48 

1.5 

0.75 

7 

1.55 

1.58 

0.99 

8 

1.77 

1.74 

1.1 

9 

1.86 

1.74 

1.18 

10 

1.82 

1.7 

1.09 

11 

1.83 

1.7 

1.15 

Table  16 

Calculated  Entropy  for  Each  Station,  by  Years  of  Data  (40  Intervals  for  Fisk, 

52  Intervals  for  Wilhelmina  and  76  Intervals  for  Clark  Comer) 

Number  of  Years  of  Data 

FISK,  MO 

WILHELMINA 

CLARK  CORNER 

1 

0.6 

1.3 

0.5 

2 

1.04 

1.89 

0.84 

3 

1.11 

1.79 

0.72 

4 

1.17 

1.75 

0.56 

5 

1.32 

1.81 

0.64 

6 

^  1.48 

1.96 

0.94 

7 

1.55 

2.04 

1.18 

8 

1.77 

2.2 

1.34 

9 

2.2 

1.45 

10 

1.82 

2.14 

1.36 

11 

1.83 

2.15 

1.45 

The  values  in  Tables  15  and  16  are  plotted  in  Figures  57  and  58,  respectively.  It  should 
be  noted  that  while  the  addition  of  more  intervals  changed  the  absolute  values  of  the 
entropy,  the  relative  relationships  were  the  same.  The  absolute  values  have  no  meaning 
except  as  compared  to  themselves,  since  the  number  of  intervals  can  affect  the  absolute 
values.  It  is  the  pattern  of  increase  that  is  important,  and  not  the  values  themselves. 
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ENTROPY 


Figure  57.  Years  of  data  vs.  incremental  entropy  (40  intervals,  Fisk,  30  intervals, 
Wilhelmina,  45  intervals,  Clark  Comer,  used  for  entropy  computations) 


Figure  58.  Years  of  data  vs.  incremental  entropy  (40  intervals,  Fisk,  52  intervals, 
Wilhelmina,  76  intervals,  Clark  Comer,  used  for  entropy  computations) 
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It  can  be  observed  in  both  Figure  57  and  Figiire  58  that  the  values  of  the  entropy 
at  the  three  stations  varies  with  time.  However,  after  eight  years  at  Wilhelmina  and  nine 
years  at  the  other  two  stations,  the  value  of  the  entropy  reaches  a  maximum.  This 
indicates  that  the  additional  data  beyond  this  point  of  maximum  entropy  is  not  decreasing 
the  bias  associated  with  the  random  variable,  sediment  transport,  nor  is  any  additional 
information  being  provided.  It  has  been  stated  that  entropy  can  be  employed  as  a 
measure  of  the  degree  of  ignorance  as  to  the  true  state  of  a  system  (Singh  and  Fiorentino, 
1992).  Given  this,  the  additional  sediment  data  beyond  year  eight  (or  nine)  was  again, 
not  adding  to  the  definition  of  the  true  values  of  the  state  variable,  computed  sediment 
transport. 

There  are  also  two  points  worth  noting  as  pertain  to  the  entropy  change  over  time. 
The  drops  in  entropy  noted  in  years  three  and  four  at  Wilhelmina  and  Clark  Comer 
indicate  that  there  is  no  additional  data  being  transferred  by  the  data  during  this  period  at 
these  stations.  In  addition,  the  drops  in  entropy  values  are  likely  an  artifact  of  the 
discretized  distribution  function,  which  can  only  approximate  the  continuous  fimction. 
The  same  conclusions  may  be  drawn  from  the  slight  fluctuations  in  entropy  values  after 
the  maximum  is  reached.  These  are  again  reflections  of  the  descritized  data.  The 
information  contained  in  the  data  set  should  and  cannot,  in  reality,  drop  with  the  addition 
of  additional  data.  It  can  and  will,  however,  reach  a  point  where  additional  data  supplies 
no  more  information.  These  are  the  maximum  values  noted  above  for  each  station. 
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Sensitivity  Analysis 


Having  developed  a  procedure  to  determine  the  optimum  length  of  data  collection 
in  the  St.  Francis  basin  based  on  historical  conditions,  a  sensitivity  analysis  was 
performed  to  evaluate  the  order  in  which  the  data  were  presented.  To  do  this,  the  entropy 
analysis  was  repeated  for  two  separate  conditions.  In  one  case,  the  order  in  which  the 
data  were  introduced  was  reversed.  That  is,  the  historic  year  one  was  introduced  as  data 
segment  1 1.  Historic  year  1 1  was  introduced  as  data  segment  one,  and  so  forth.  This  was 
performed  for  all  three  stations  and  the  results  (maximums  boldface)  are  shown  in 
Table  17. 


Table  17 

Calculated  Entropy  for  Each  Station,  by  Years  of  Data,  with  the  Historic 

Order  of  Occurrence  Reversed  (40  Intervals  for  Fisk,  52  Intervals  for 
Wilhelmina  and  76  Intervals  for  Clark  Comer) 

Number  of  Years  of  Data 

FISK,  MO 

WILHELMINA 

CURK  CORNER 

1 

1.21 

1.32 

1.07 

2 

1.36 

1,3 

0.91 

3 

1.77 

1,63 

1.36 

4 

2.09 

2.07 

1.59 

5 

2.07 

2.14 

1.69 

6 

2.1 

2.23 

1.72 

7 

2.1 

2.18 

r  1.69 

8 

2 

2.17 

1.54 

9 

1-91 

2.1 

1.41 

10 

1.9 

2.19 

1.45 

11 

1.83 

2.15 

1.45 
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It  can  be  seen  that  reversing  the  historical  order  of  occurrence  of  the  data  changed 
the  point  of  maximum  entropy.  In  each  case,  the  maximum  occurred  after  six  years  of 
data  were  introduced.  While  different  from  the  results  of  the  historical  data  analysis,  the 
results  were  consistent.  Based  on  this  portion  of  the  sensitivity  analysis,  it  appears  that 
the  approach  is  still  valid.  One  would  expect  different  results  with  different  data  sets. 

The  results  are,  however,  consistent  within  the  basin.  They  reflect  a  real  pattern  of  events 
(if  reversed)  that  might  occur.  Hydrologic  conditions  are  marked  by  cycles  of  drought, 
flood  and  “normal”  conditions.  Therefore,  this  reversed  data  still  represent  a  possible 
sequence  of  events.  The  procedure  would  still  determine  the  optimum  length  of  record  to 
present  the  maximum  amount  of  information,  in  this  case  six  years. 

Finally,  the  data  were  introduced  in  a  totally  random  fashion.  The  data  sets 
generated  in  this  manner  bore  no  resemblance  to  one  another  from  one  station  to  the  next. 
In  this  case,  the  resulting  sediment  transport  was  not  generated  by  real  global  and 
regional  weather  patterns  occurring  over  the  basin,  but  rather  as  disjointed  and  random 
events.  Results  derived  from  this  action  would  not  be  expected  to  be  uniform  for  the 
three  stations.  Results  are  shown  in  Table  18. 

In  this  hypothetical  world,  the  results  are  different  from  the  previous  two  cases.  It 
cannot  be  determined  if  the  values  at  Fisk  and  Clark  Comer  have  reached  their  maximum 
value.  It  appears  from  the  rate  of  increase  in  their  entropy  values  that  these  station  are  at 
or  near  the  maximum,  but  more  data  would  be  necessary  to  confirm  this.  At  Wilhelmina, 
the  maximum  was  reached  after  the  sixth  data  segment  was  introduced.  Data  introduced 
in  this  fashion  is  of  interest  only  from  a  theoretical  standpoint.  In  reality,  the  data  at  the 
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Table  18 

Calculated  Entropy  for  Each  Station,  by  Years  of  Data,  with  the  Historic 
Order  of  Occurrence  Introduced  Randomly  for  Each  Station.  (40  Intervals  for 
Fisk,  52  Intervals  for  Wilhelmina  and  76  Intervals  for  Clark  Corner) 


Number  of  Years  of  Data 

FISK,  MO 

WILHELMINA 

CLARK  CORNER 

1 

0.9 

0.9 

0.64 

2 

0.93 

1.28 

0.52 

3 

1.09 

1.51 

0.85 

4 

1.26 

1.8 

0.78 

5 

1.22 

1.9 

0.82 

6 

1.33 

r  2.16 

1.13 

7 

1.6 

2.12 

0.99 

8 

1.65 

2.13 

1.18 

9 

1.7 

2.06 

1.31 

10 

1.79 

2.12 

1.42 

11 

1.83 

2.15 

1.45 

three  station  would  be  somewhat  inter-dependent  in  that  the  hydrology,  temperatures  and 
weather  conditions  would  be  similar  for  each  station  in  a  given  year.  This  random 
approach  nullifies  that  relationship.  In  spite  of  this,  there  are  no  spurious  results  that 
would  invalidate  the  procedure. 

Having  determined  that  the  procedure  can  indeed  identify  the  length  of  record 
necessary  to  provide  the  maximum  amoimt  of  information  at  the  three  stations,  it  next 
seemed  prudent  to  determine  the  impacts  of  collecting  less  than  the  optimal  amount  of 
data.  For  example,  suppose  that,  based  on  the  research  herein  using  the  historical  data, 
one  might  believe  that  eight  or  nine  years  or  data  were  necessary  at  any  point  in  the 
St.  Francis  basin.  However,  suppose  that  due  to  time  and/or  economic  constraints,  that 
one  had  only  five  years  available  to  collect  data  at  another  point  within  the  basin.  What 
would  be  the  impact?  To  answer  this  question,  the  maximum  entropy  values  were  used 
as  a  base  value.  Next,  the  incremental  entropy  was  divided  by  this  base  value.  This 
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yields  the  amoxint  of  the  total  information  present  at  the  end  of  each  data  segment.  These 
values  are  presented  as  a  percent  in  Table  19  for  the  40  segment  Fisk  analysis,  the 
52  segment  Wilhelmina  analysis  and  the  76  segment  Clark  Comer  Analysis.  Since 
information  cannot  increase  beyond  the  point  of  maximum  entropy,  calculations  beyond 
this  point  would  be  meaningless  and  these  are  not  provided.  These  data  are  shown 
graphically  in  Figure  59. 


Table  19 

Quantity  of  Information  Supplied  by  Station  with  Increasing  Size  of  the  Data  Set, 
as  a  Percent  of  the  Total 

Years  of  Data 

Fisk 

Wilhelmina 

Clark  Comer 

1 

32 

59 

34 

2 

56 

85 

58 

3 

60 

85* 

58* 

4 

63 

85* 

58* 

5 

71 

85* 

58* 

6 

80 

89 

64 

7 

83 

93 

81 

8 

95 

100 

92 

9 

100 

- 

100 

10 

- 

- 

- 

11 

- 

*  Note,  actual  amount  of  information  cannot  decrease,  as  calculations  would  seem  to  indicate. 
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Figure  59.  Percent  of  information 


The  results  of  this  calculation  show  that  if  only  five  years  of  data  could  be 
collected,  these  data  would  provide  between  58  and  71  percent  of  the  maYimnm 
information  one  would  get  by  collecting  the  recommended  amount  of  data.  Using  these 
data  and  the  regression  techniques  above,  one  would  select  the  Yang-HEC6  transport 
function  at  Fisk,  the  same  decision  reached  with  the  full  information  as  shown  in 
Table  10.  At  Wilhelmina,  the  Brownlie  transport  function  would  be  selected  based  on  the 
five  years  of  data,  also  the  same  decision  reached  with  the  full  information  (see  Table  12). 
However,  at  Clark  Comer,  one  would  select  the  Ackers- White  HEC6  transport  fimction, 
and  with  some  confidence.  It  is,  after  all,  ranked  number  one  at  the  end  of  years  2-5. 
However,  once  one  has  the  full  information,  the  Ackers- White  function  slips  to  10th  and 
the  best  selection  is  the  Toffaleti-Schoklitsch  transport  function,  which  was  ranked  3rd 
after  year  five  (see  Table  14).  It  should  also  be  noted  that  after  year  five,  71  and 
85  percent  of  the  maximum  information  is  available  at  Fisk  and  Wilhehnina,  while  only 
58  percent  is  available  at  Clark  Comer.  Therefore,  at  the  stations  with  a  high  degree  of 
information  available,  the  correct  decision  would  be  made.  At  Clark  Comer,  with  barely 
half  of  the  information  available,  an  abysmal  decision  would  be  made  that  looked  good 
based  on  other  commonly  used  selection  criteria.  This  analysis  would  suggest  that  for 
new  stations  near  Fisk  and  Wilhelmina,  5  years  of  record  would  be  sufficient,  though  not 
optimum.  New  stations  in  the  vicinity  of  Clark  Comer  require  at  least  eight  years  of 
record  to  insure  enough  information  is  available  to  make  the  correct  decision  as  to  which 
transport  function  to  use.  In  any  case,  the  new  station  data  should  be  checked  with  the 
procedure  presented  in  the  next  chapter  to  determine  the  proper  length  of  the  data 
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collection  effort. 


5  A  Procedure  for  Evaluating  the 
Adequacy  of  Sediment  Data  Sets 


The  Procedure 


The  following  procedure  is  to  be  used  for  existing  sediment  data  collection 
stations.  The  initial  step  is  to  select  candidate  sediment  transport  functions  for  use  in  the 
analysis.  In  this  dissertation,  ten  separate  functions  were  used  for  research  purposes.  In 
practical  application,  two  or  three  can  be  selected  based  on  their  applicability  to  the 
stream  to  be  analyzed.  These  selection  criteria  for  the  transport  functions  are  available 
from  several  sources  (Vanoni,  1975;  Thomas  et  al.,  1992).  Once  this  selection  has  been 
made,  the  data  may  require  augmentation  by  computing  water  surface  slopes  or 
compiling  grain  size  distribution  data  information.  The  sediment  computations  can  then 
be  made.  These  computations  will  compute  sediment  transport  for  each  observation  based 
on  hydraulic  and  sediment  parameters  measured  during  the  observation.  Computer 
packages  such  as  the  U.  S.  Army  Corps  of  Engineers  SAM  (Thomas  et  al.,  1992) 
software  greatly  ease  these  computations. 
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The  sediment  computations  should  be  conduced  in  segments.  The  first  segment 
will  contain  only  one  year  of  hydraulic  and  sediment  data  with  which  to  make 
calculations.  Each  successive  segment  will  add  an  additional  year  of  data,  so  that 
segment  two  includes  data  from  years  one  and  two,  segment  three  data  from  years  one, 
two  and  three  and  so  forth.  A  regression  analysis  may  then  be  performed  between  the 
computed  and  observed  sediment  discharge.  The  transport  function’s  performance  can  be 
ranked  for  each  segment  by  comparing  the  coefficients  of  determination.  The  rank  of 
these  coefficients  can  then  be  plotted  versus  time.  The  point  at  which  the  values  of 
relative  rank  cease  to  change  will  identify  an  initial  range  of  minimum  data  necessary  to 
allow  meaningful  sediment  transport  computations  to  be  made.  Should  no  pattern 
develop  and  the  rank  continue  to  change  imtil  all  data  are  added,  then  sufficient  data  are 
not  available  to  define  the  sediment  transport  relationship  at  that  station,  or  different 
transport  functions  should  be  selected  and  the  process  repeated. 

Next,  the  computed  values  should  be  fitted  to  a  distribution  function  and  the  total 
entropy  calculated  for  the  data  with  succeeding  years  of  data  incrementally  introduced.  If 
the  entropy  ceases  to  increase  at  some  point  and  remains  the  same  thereafter,  then  that 
point  defines  the  amoimt  of  data  necessary  to  define  the  sediment  transport  relationship. 
Slight  fluctuations  in  the  value  after  this  point  are  to  be  expected  as  artifacts  of  the 
discretized  distribution  function.  If  no  such  point  is  reached,  then  insufficient  data  have 
been  collected  and  the  data  collection  program  should  continue.  If  this  point  lies  beyond 
the  point  suggested  by  the  regression  analysis,  then  the  entropy  driven  point  should  take 
precedence. 
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Once  such  a  point  has  been  reached,  the  percentage  of  information  present  after 
addition  of  each  segment  of  data  can  be  computed  based  on  a  percentage  of  the  maximum 
entropy  value.  From  this,  inferences  can  be  made  for  other  locations  within  the  basin 
where,  perhaps,  less  data  are  available  as  to  what  percent  of  the  information  is  provided 
by  a  data  set  of  less  than  optimum  length. 

It  occurs  to  the  author  that  there  are  conditions  under  which  one  may  wish  to 
continue  data  collection  even  if  the  procedure  outlined  above  indicates  it  is  no  longer 
necessary.  If  conditions  are  expected  or  could  occur  which  would  significantly  change 
the  distribution  of  the  data,  then  one  would  want  to  continue  the  data  collection.  This 
would  be  necessary  because  a  significant  change  in  the  distribution  could  increase  the 
entropy.  One  such  instance  would  be  if  significant  widespread  changes  in  land  use  are 
anticipated.  This  could  affect  the  sediment  water-sediment  discharge  relationships  and 
alter  the  distribution  of  previously  collected  data.  Additional  data  would  add  information 
about  the  sediment  relationships,  in  this  case.  Other  likely  causes  would  be  significant 
water  diversions,  dam  construction  or  removal,  or  natural  disasters  such  as  volcanic 
activity. 

The  research  described  herein  was  conducted  at  several  sites  within  the  St.  Francis 
River  basin  that  were  selected  to  reflect  rather  different  conditions.  It  was  believed  that 
these  varied  sufficiently  so  as  to  essentially  be  entirely  different  in  nature.  The  Fisk,  Mo., 
Station  is  located  just  below  Wappapello  dam  in  a  rather  stable  section  of  the  river,  with 
an  average  drop  in  channel  grade  of  one  foot  per  mile,  higher  than  that  found  in  the  lower 
reaches.  The  bed  material  samples  also  varied  significantly  at  this  station,  resulting  in  the 


145 


very  low  coefficients  of  determination  as  shown  in  Table  6.  The  Wilhelmina  station  is 
located  in  a  highly  unstable  man-made  reach  of  the  river  which  requires  channel  clean  out 
periodically  to  restore  the  hydraulic  conductivity.  Slopes  are  relatively  high  here, 
averaging  5.08  feet  of  drop  per  mile  of  chaimel,  as  the  cutoff  reduced  the  natural  channel 
length  by  7  miles.  This  reach  experiences  continuous  shoaling  and  displays  a  high 
sediment  transport  load.  Even  though  this  was  a  very  active  reach  for  sediment  transport, 
the  coefficients  of  determination  were  higher  than  those  computed  for  Fisk,  as  shown  in 
Table  7.  The  Clark  Comer  station  is  located  in  another  man-made  reach  but  is  far  more 
stable.  This  reach  has  relatively  mild  slopes,  averaging  0.63  feet  per  mile,  and  is  just 
above  the  section  of  the  river  that  experiences  backwater  effects  from  high  Mississippi 
river  stages  (USACE,  1992).  Thus,  the  procedure  described  above  was  developed  under 
rather  different  regimes  and  should  be  robust  enough  to  account  for  changing  conditions 
within  parameters  normally  observed  in  most  agricultural  basins. 


6  Summary  and  Conclusions 


One  of  the  main  purposes  for  collecting  sediment  data  is  to  provide  sufficient 
data  to  select  a  transport  function  that  best  represents  the  observed  data.  This  allows  one 
to  make  predictions  as  to  future  conditions,  analyze  flows  that  are  different  from  or 
outside  of  the  observed  data  set,  and  evaluate  system  responses  to  proposed  changes  in 
the  system.  These  changes  may  be  man  induced  such  as  land  use  changes,  structures 
placed  in  the  channel  or  other  forms  of  channel  modification.  Tfre  changes  may  also  be 
due  to  natural  causes,  though  unless  the  natural  occurrence  is  of  a  violent  nature,  such  as 
an  earthquake,  volcanic  action  or  forest  fire,  the  time  scale  is  usually  too  long  to  be  of 
significance  for  engineering  purposes.  The  procedure  presented  in  the  previous  chapter 
will  allow  the  sediment  engineer  to  assess  existing  data  sets  and  determine  if  they 
adequately  describe  the  sediment  transport  characteristics  of  the  stream.  If  not,  inferences 
can  be  made  based  on  the  degree  of  change  in  the  ranking  of  the  transport  functions  and 
the  rate  of  change  of  the  entropy  as  to  how  well  the  data  describe  the  system. 

It  is  common  practice  to  compare  observed  sediment  data  with  several  transport 
functions  and  use  the  one  that  best  matches  the  data  (Julien,  1995).  Various  methods 


have  been  used  in  the  past.  Visual  comparison  of  the  computed  versus  the  measured 
sediment  discharge  was  used  by  Vanoni  et  al.  (1961).  This  can  be  demonstrated  by 
looking  at  the  plots  in  Figures  6-11.  The  drawbacks  to  this  approach  are  obvious.  While 
the  transport  functions  that  match  neither  the  trend  or  the  vales  are  obvious,  it  is  difficult 
to  distinguish  between  those  that  cluster  near  the  observed  data.  Other  researchers 
(White  et  al.,  1973;  Yang  and  Wan,  1991)  have  applied  the  discrepancy  ratio  (the  ratio  of 
the  observed  to  the  calculated  values)  as  a  measure  of  which  transport  function  best 
matches  the  observed  data.  Since  the  ratios  are  almost  never  “one”,  it  is  difficult  to 
interpret  the  differences,  which  are  numbers  such  as  3.6  or  0.45,  at  a  glance.  A  different 
approach  has  been  proposed  by  Gomez  and  Church  (1989).  This  method  uses  the 
analysis  of  variance  (ANOVA)  method  to  compare  local  bias  and  mean  bias  values 
between  the  computed  and  observed  sediment  transport.  However,  this  method  also 
displays  results  as  a  distribution  of  the  ratio  of  the  computed  values  to  those  of  the 
observed.  Again,  the  results  require  some  detailed  analysis  to  interpret. 

The  regression  method  used  in  this  dissertation  to  rank  the  sediment  transport 
functions  and  select  the  best  one  is  unique.  While  regression  methods  have  been  applied 
previously  by  other  researchers,  they  have  never  been  applied  incrementally  to  segments 
of  the  data  set.  This  method  requires  that  a  transport  function  be  the  best  ranked 
incrementally  over  the  data  set.  Furthermore,  this  method  specifies  that  a  transport 
function  also  be  the  best  performer  coincidentally  with  the  maviTTmm  amount  of 
information.  The  incremental  analysis  points  out  the  dangers  involved  with  selecting  a 
transport  function  based  on  a  finite  set  of  data,  especially  if  the  data  set  is  not  providing 
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the  maxinnim  amount  of  information.  The  results  of  this  procedure  are  easy  to  interpret, 
either  from  tables  or  graphically. 

The  regression  method  applied  herein  was  first  performed  with  the  full  data  set  for 
each  station.  Based  on  the  analyses  at  the  three  stations,  a  different  transport  function 
best  matched  the  observed  data  at  each  station.  Furthermore,  each  of  the  number  one 
ranked  equations  had  a  different  theoretical  basis.  At  Fisk,  MO,  the  Yang-HEC6 
transport  function  was  the  highest  ranked.  It  is  based  on  stream  power  and  has  been 
modified  by  the  WES  to  calculate  sediment  transport  by  size  class.  At  Wilhelmina,  Mo., 
the  Brownlie  Function  best  matched  the  observed  data.  This  function  is  an  empirical 
formulation  based  on  regression  analysis.  At  Clark  Comer,  AR,  the  hybrid  Toffaleti- 
Schoklitsch  function  best  matched  the  observed  data.  This  formula  calculates  the 
sediment  transport  using  both  the  methods  of  Toffaleti  and  those  of  Schoklitsch  and  uses 
the  higher  value.  This  hybrid  function  only  marginally  outperformed  the  basic  Toffaleti 
function  at  this  station,  as  indicated  by  comparing  the  coefficients  of  determination.  The 
value  based  on  the  full  data  set  for  the  hybrid  was  78.46  versus  78.44  for  the  basic 
function.  This  indicates  little  was  gained  by  adding  the  bed  load  computations.  This 
function  is  based  on  the  advection-diffusion  approach  of  Einstein.  The  Schoklitsch 
function  is  an  empirical  formula  that  is  primarily  recommended  for  bed  load 
computations.  Having  noted  the  variation  of  functions  and  their  theoretical  bases,  overall 
the  Brownlie  function  was  the  best  performer.  It  ranked  3, 1  and  3  at  Fisk,  Wilhelmina 
and  Clark  Comer  respectively.  The  Toffaleti-Schoklistch  function  ranked  2,  6  and  1  and 
the  Yang-HEC6  function  ranked  1, 4  and  4  respectively.  Brownlie’s  average  finish  was 
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2.33  compared  to  3.0  for  both  Yang-HEC6  and  Toffaleti-Schoklistch.  All  three  of  these 
were  excellent  performers.  This  was  not  unexpected  given  that  they  quite  often  finish 
very  high  when  compared  to  other  transport  functions  using  actual  river  data  (Yang  and 
Molinas,  1991;  Brownlie,  1981). 

It  has  been  shown  that  reliance  totally  on  regression  methods,  if  the  data  record  is 
of  insufficient  length,  may  lead  to  selection  of  the  incorrect  transport  function  and/or 
validate  too  short  a  sampling  period.  For  example,  based  solely  on  regression  analysis, 
for  the  stations  analyzed  in  this  report,  one  might  look  at  one  year  of  data  and  select  the 
Colby-HEC6  transport  function  for  Fisk,  Mo.  However,  based  on  the  full  eleven  years  of 
record,  the  Colby-HEC6  function  finished  seventh  out  of  ten  at  Fisk.  If  one  looked  at 
four  years  of  record  at  Fisk,  one  would  select  the  Yang-HEC6  transport  function.  Indeed, 
this  is  the  number  one  rated  function  even  after  eleven  years  of  data  have  been  collected. 
However,  looking  at  the  entropy  analysis  reveals  that  four  years  provide  only  63  percent 
of  the  mformation  available  after  nine  years  at  Fisk.  Use  of  relative  rank  of  the  transport 
functions  as  determined  by  regression  analysis  is  at  best  good  for  determining  a  broad 
range  of  values  within  which  the  true  minimum  record  length  resides.  The  regression 
methods  presented  herein  are  good,  however,  for  determining  the  transport  function  that 
best  fits  the  observed  data.  Simply  looking  at  the  base  data  from  which  the  functions 
were  developed  can  lead  to  very  erroneous  selections.  For  instance,  based  on  the  average 
size  of  the  bed  sediment  and  sand  bed  nature  of  the  St.  Francis  River,  the  Ackers- White 
transport  function  would  be  expected  to  be  a  good  predictor.  To  the  contrary,  the  single 
grain  size  function  ranked  overall  10,  9,  and  10  at  Fisk,  Wilhelmina  and  Clark  Comer 
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respectively.  The  multiple  grain  size  version  did  only  marginally  better  with  ranks  of 
4, 10,  and  9  overall.  As  discussed  earlier,  this  was  imdoubtedly  due  to  variations  in  the 
bed  material.  The  Colby-HEC6  function  was  also  a  poor  performer.  This  function 
ranked  7, 2,  and  8,  respectively,  at  Fisk,  Wilhelmina  and  Clark  Comer.  This  function  was 
developed  from  flume  data  and  most  values  in  the  10  foot  depth  and  all  in  the  100  foot 
depth  range  were  extrapolated  from  flume  data.  Therefore,  this  function  may  not  be  the 
best  for  larger  natural  streams  that  may  often  have  flow  depths  greater  than  10  feet, 
including  the  St.  Francis  River.  The  incremental  regression  method  shows  which 
functions  are  the  best  performers  over  a  broad  range  of  data. 

The  change  in  order  of  ranking  of  the  transport  equations  as  data  are  added  seems 
to  have  some  significance.  Using  incremental  analysis  of  the  data,  if  the  order  of  ranking 
of  the  transport  functions  are  similar,  Avith  only  a  small  change  in  order  toward  the  end  of 
the  data  set,  then  the  data  woAild  seem  to  be  nearly  adequate.  If,  on  top  of  these 
conditions,  the  rate  of  change  in  entropy  is  very  small  or  zero,  then  additional  proof  is 
offered  that  the  data  set  is  adequate.  The  critical  conditions  for  selection  of  the  best 
transport  equation  are; 

1.  The  range  of  key  variables,  such  as  velocity,  depth  and  bed  sediment  size,  are 
within  the  range  recommended  for  the  transport  function. 

2.  The  top  ranked  function  remains  so  over  several  data  segments. 

3.  Relatively  stable  relative  ranking  of  the  best  performers  over  a  portion  of  the 
data  set,  particularly  as  the  optimum  length  of  the  data  set  is  approached. 
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4.  The  optimal  length  of  data  collection  is  determined  using  entropy  principles 
and  the  selected  transport  function  should  be  top  ranked  at  this  point. 

If,  on  the  other  hand,  the  rank  order  of  the  equations  is  varying  wildly  or  the  rate  of 
change  in  the  entropy  is  not  zero,  then  strong  evidence  exists  that  the  data  set  is  not 
sufficient  to  describe  the  sediment  transport  characteristics  of  the  stream. 

The  principle  of  entropy  has  been  successfully  applied  to  determine  the  optimum 
length  of  sediment  data  collection  at  three  stations  in  the  St.  Francis  river  basin.  This 
optimum  length  was  nine  years  at  two  stations  (Fisk,  Clark  Comer)  and  eight  years  at 
Wilhelmina,  based  on  the  historical  data.  This  is  the  first  definitive  guidance  developed 
as  to  how  much  data  need  be  collected  to  provide  the  most  information  on  sediment 
transport  at  a  given  location.  With  the  historic  order  of  data  reversed,  the  optimum  length 
of  data  collection  was  6  years  and  was  consistent  at  each  station.  Analysis  of  this  showed 
the  difference  in  results  were  due  to  the  order  in  which  the  distribution  intervals  were 
filled.  After  nine  years  at  Fisk,  there  were  85  observations  placed  in  17  of  40  intervals. 
Forty  five  of  these  were  clustered  in  the  first  interval.  After  6  years  at  Fisk,  with  the 
order  of  data  reversed,  there  were  47  observations  placed  in  13  of  40  intervals.  Eighteen 
of  these  were  in  the  first  interval.  Therefore,  though  the  historic  record  had  38  more 
observations  than  the  reverse  order  data  set,  27  of  these  were  in  the  first  interval  and 
added  little  information.  A  similar  pattern  was  observed  at  Clark  Comer,  in  the  lower 
end  of  the  basin.  After  nine  years,  there  were  68  observations  placed  in  13  of 
76  intervals.  Forty  two  of  these  were  clustered  in  the  first  interval.  After  6  years  at  Clark 
Comer,  with  the  order  of  data  reversed,  there  were  43  observations  placed  in  12  of 
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76  intervals.  Nineteen  of  these  were  in  the  first  interval.  Therefore,  though  the  historic 
record  had  25  more  observations  than  the  reverse  order  data  set,  23  of  these  were  in  the 
first  interval  and  added  little  information.  Therefore,  the  procedure  is  consistent  and 
awards  short  data  sets  Avith  better  distribution  over  larger,  more  poorly  distributed  data 
sets. 

Therefore,  the  basic  question  that  has  been  answered  and,  indeed,  is  the  topic  of 
this  research,  is  how  much  data  are  needed?  Historically,  the  answer  has  been  “as  much 
data  as  can  be  afforded  within  study  budget  constraints  of  time  and  money”.  Clearly  if 
this  method  has  ever  provided  adequate  data  it  was  serendipity  as  in  most  all  cases  either 
too  little  or  too  much  data  would  be  collected.  The  above  outlined  procedure  (Chapter  5) 
allows  the  sediment  engineer  to  evaluate  existing  data  sets  as  to  their  adequacy,  that  is, 
answer  the  question  “Do  I  have  enough  data?”.  It  also  allows  one  to  select  the  best 
sediment  transport  equation  given  the  available  data,  and  make  a  determination  as  to  how 
well  the  equation  will  describe  the  sediment  transport  within  the  stream  reach. 

Predictions  and/or  evaluations  can  then  be  made  vvdth  some  degree  of  surety  as  opposed 
to  assuming  data  available  are  adequate  and  that  the  selected  sediment  transport  equation 
is  the  best  available.  The  entropy  procedure  also  provides  a  quantitative  indication  as  to 
the  amount  of  information  available  for  the  predictive  system.  The  procedure  has  been 
proven  robust  through  sensitivity  analyses.  While  randomly  introducing  the  data 
destroyed  the  pattern  observed  between  historical  data  collection  stations  and  the  same 
data  introduced  in  reverse  order,  it  none-the-less  was  still  useful  in  determining  if  the 
maximum  amount  of  information  was  present.  It  has  also  been  proven  that  the 
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distribution  of  the  data  is  more  important  than  the  qiiantity  of  data. 

Thus,  one  can  use  the  methods  of  this  report  and  determine  if  sufficient  data  are 
available.  If  not,  then  some  notion  as  to  what  percent  of  the  maximum  information  is 
available  can  be  made.  This  can  be  determined  by  looking  at  the  rate  of  increase  in  the 
entropy.  If  it  is  increasing  rapidly,  then  the  rate  of  information  addition  is  high  and  one 
additional  year  of  data  would  be  more  valuable.  If  the  rate  of  change  is  low,  then 
additional  data  will  likely  not  add  much  additional  information. 

The  procedure  developed  herein  has  application  in  real  world  decision  making. 
For  example,  the  original  St.  Francis  river  sediment  data  collection  program  consisted  of 
24  stations.  This  number  was  reduced  in  1983  and  again  in  1985.  These  reductions  were 
dictated  by  economic  necessity.  In  actuality,  the  decisions  were  made  based  on  the 
perceived  goodness  of  the  data  collected  and  importance  of  each  station’s  location.  Using 
the  method  of  this  report,  a  quantitative  basis  could  have  been  established  for  removing 
or,  keeping,  individual  stations.  It  is  possible  that  in  1985  a  number  of  stations  could 
have  been  dropped  based  on  the  findings  at  in  this  report.  A  case  could  also  have  been 
made  to  continue  sampling,  checking  yearly,  until  the  optimum  length  had  been  reached 
at  individual  stations.  This  could  be  done  in  the  following  manner.  If  one  were  to 
conduct  an  incremental  analysis  of  information,  normalized  to  a  base  year  taken  to  be  the 
proceeding  year,  then  the  rate  and  relative  value  of  information  added  by  taking 
additional  data  could  be  determined.  For  example,  let  us  look  at  the  Wilhelmina  station. 
Using  the  values  jfrom  Table  16.  The  incremental  entropy,  or  information  added,  can  be 
calculated  as 
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1.89-1.30 

- =.45 

1.30 

Continuing  this  process  yields  the  following  results,  as  shown  in  Table  20.  A  case  could 
be  made  to  continue  collecting  data  prior  to  year  9.  After  that  time,  cessation  of  data 
collection  at  this  station  could  be  justified.  Of  course,  one  may  not  know  the  maYimnm 
value  has  been  reached  until  the  optimum  number  of  years  plus  one  of  data  have  been 
collected.  This  is  because  to  identify  the  maximum  value  would  require  at  least  one  more 
year  of  data,  as  can  be  inferred  from  Table  16. 


Table  20 

Incremental  Entropy  for  Each  additional  Year  of  Data,  Wilhelmia,  Mo. 

Incremental  Entropy  From _ year  To _ year 

Incremental  Entropy 

1-2 

0.45 

2-3 

0.00* 

3-4 

0.00* 

4-5 

0.03 

5-6 

.08 

6-7 

.04 

7-8 

.08 

8-9 

0.00** 

*  Although  the  calculated  values  is  less  than  zero,  Information  cannot  be  deleted  by  additional  sampling. 
These  values  are  due  to  the  discretized  distribution. 

**  Actual  maximum  entropy  location 
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The  procedure  presented  herein  could  also  be  applied  to  assist  in  designing  a  dgta 
collection  program.  If  data  were  available  in  a  basin  and  had  previously  been  analyzed 
using  the  methods  of  this  dissertation,  then  one  could  make  an  educated  estimate  of  the 
length  of  data  needed  at  a  new  location  within  the  basin.  The  data  collected  at  the  new 
location  could  be  checked  with  the  incremental  method  outlined  above  and  some  notion 
provided  as  to  the  amount  of  information  gained  with  each  additional  year  of  data 
collected. 

Of  course,  like  any  new  procedure,  this  one  needs  to  be  applied  by  other 
researchers  to  streams  other  than  the  St.  Francis  River.  It  has  been  rather  rigorously 
applied  in  the  St  Francis  basin  and  can  be  used  to  evaluate  data  in  this  basin.  It  is 
recommended  that  this  procedure  be  applied  to  a  number  of  streams.  In  so  doing,  it  is  my 
belief  that  a  body  of  knowledge  can  be  assembled  that  will  offer  insight  as  to  just  how 
much  data  are  required  for  a  given  type  of  stream,  be  it  a  medium  sized  alluvial  basin 
such  as  the  St.  Francis  River  basin,  or  one  with  different  material  types  and  slope 
characteristics.  Over  time,  it  may  come  to  pass  that  the  length  of  data  for  each  class  of 
stream  can  be  determined  fairly  accurately,  from  climate  and  basin  characteristics  much 
as  Manning’s  “n”  values  are  now  determined  by  a  cursory  inspection  of  a  channel  reach. 
However,  it  should  be  noted  that  due  to  sensitivity  of  the  procedure  to  the  order  in  which 
the  data  are  added,  that  the  author  would  not  recommend  directly  transferring  results  from 
one  basin  to  another  until  further  research  has  been  done  on  the  performance  of  the 
procedure  on  similar  but  geographically  separated  basins. 
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In  addition  to  the  above  procedures,  the  author  has,  in  pursuit  of  this  research, 
assembled  a  truly  unique  set  of  sediment  data.  These  data,  much  of  which  are 
unpublished  elsewhere,  are  an  extremely  valuable  asset.  They  can  serve  as  a  basis  to 
further  the  research  described  herein,  or  to  provide  other  researchers  with  an  invaluable 
data  source  to  further  the  knowledge  of  sediment  transport.  Hopefully,  the  use  of  this 
data  may  refine  existing  methods  and  procedures  and/or  lead  to  new  findings  which 
further  the  field  of  Sedimentation  Engineering. 
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APPENDIX  A 


SAMPLE  COMPUTER  INPUT/OUTPUT 


SAMPLE  SAS  OUTPUT 


Chronological  list  of  data  for  Fisk,  Wilhelmina  and  Clark  Comer 
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14  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 
4  4  4  4  4 

•  4  4  4  4 

4  14  4  4 
4  4  4  4  4 

1  4  4  4  4 

4  4  1  4  4  1 

4  4  4  4  4  1 

4  4  4  4  4  3 

•  44143 

•  4  4  4  4  3 
«  4  4  4  I  3 

•  4  4  4  4  4 

•  4  4  4  4  1 
4  4  4  4  •  I 
4  I  4  4  4  4 

•  4  4  4  4  4 

•  4  4  4  4  4 

•  4  4  4  I  4 

•  4  4  4  4  4 
I  4  I  4  4  4 
1  4  4  4  4  4 

9  4  4  4  4  4 
I  •  4  4  •  4 
I  «  4  4  4  4 

10  4  4  4  4 
I  0  4  4  0  • 

3  0  4  4  4  4 
0  4  4  0  4 
4  4  4  4  4 
0  0  4  4  4 
0  4  4  4  4 

•  0  I  4  4 

•  4  4  4  4 

•  0  4  4  4 

0  4  4  4  4 
0  0  I  4  4 

0  0  4  4  1 

0  I  •  I  4 

0  0  4  4  4 

0  0  4  4  4  1 

4  0  4  4  4  1 

0  4  4  4  4  3 


9  4  4  4  4  • 

9  4  4  4  4  1 

9  4  4  4  4  4 

I  4  I  4  4  4 

14  4  4  4  4 

I  4  4  4  4  4 

•  4  4  4  4 

4  4  4  4  4 

•  I  4  4  4 

4  4  4  4  •  7 

4  4  4  4  4 

4  4  4  I  4 

4  4  4  4  4 

4  4  4  4  4 

•  4  4  4  1 

4  4  4  4  4 

4  4  4  4  4 

4  4  4  1  4 

4  4  4  4  4  i 

•  4  4  4  4  3 

•  4410  3 

•  4  4  4  4  3 

•  •  4  4  4  3 

4  4  4  4  39  4 

4  4  4  4  4  4 

4  4  4  4  1  0 

4  4  4  4  4  4 

0  I  4  4  0  4 

0  4  4  4  1  4 

9  4  I  4  0  0 

9  4  4  4  0  0 

9  4  4  4  1  0 

9  4  4  4  0  0 

)  0  4  0  0  • 

14  4  0  0  4 

•  4  4  0  4 

•  4  4  4  0 

4  4  4  0  0 

4  4  4  4  0 

•  I  4  0  4 

4  4  4  0  4 

4  4  4  0  0 

4  4  4  0  4 

0  0  0  0  1 

1  0  4  0  0 

0  0  0  0  0 

•  4  4  I  • 

4  4  4  0  4 

4  4  4  4  0 

4  4  •  4  4 


9N  11  n  1  1 
•  m  11  21  1 

wo  M  «  M  2 

N  II  11  t  1 
N  10  44  1  2 
4  m  1$  4  4 

N  11  M  4  4 

90  11  10  S  S 

K  11  IS  2  1 

91  14  47  0  S 

0  IN  n  0  0 

4  IN  11  4  1 

M  tt  14  S  S 
«  11  6  2  I 

0  «  44  1  1 
0  11  IS  S  2 
1  11  6  1  I 
9  If  6  2  2 
9  11  14  1  I 

9  IN  e  S  1 
9  11  M  S  4 
I  IN  M  t  1 
I  «  e  1  1  3 


I  11  47  1  1 

9  11  N  1  1 

3  11  00  1  I 

3  11  77  I  1 

tt  St  t  1 
11  U  S  2 
It  IS  S  < 
11  C  1  1 
tt  07  1  1 
tt  $4  4  4 
17  6  11  10 
n  tt  t  4 


9  11  O  2  1 

9  11  40  H  2 

W  R  2  2 

IN  tt  7S  If 
11  12  IS  S 

11  N  U  S 

e  74  S  4 

17  71  $  4 

11  U  2  1 

tt  01  1  1 


1  Ut  OCS  1.71 
1  Ut  an  i.n 
i  UT  tu  1.14 
1  Ul  120  1.H 
S  le  ISM  2.6 
4  211  1721  2.01 

4  MO  Ttt  1.10 

s  te  iNo  i.n 

1  67  SON  2.70 

5  IN  CS  1.01 
•  170  OM  i.n 
1  241  UN  2.11 
S  2S1  2120  2.17 
I  220  UN  2.11 
I  2S  MM  2.6 

t  2S  tlM  I.e  1 

1  m  se  i.e 

2  Ul  221  1.6 

4  Ml  SIN  I.S1  3 
1  140  UM  2.07  1 
4  2M  SON  2.M  1 
1  IC  UM  2.U  1 

1  241  UM  2.U  1 

2  2St  SIM  2.R  t 
1  UO  ISTf  2.1S  t 
1  2M  MSI  2.S0  11 
1  2U  t$71  2.6  li 
I  2SS  UM  2.11  1^ 
1  2U  lUO  2.N  1 
t  U2  in  2,27  1 
I  U7  ISM  2.6  12 
I  241  2041  t.N  11 

3  Ul  1111  1.77  11 
>  2M  tIN  2.6  1 

171  Ul  0.71  4 
Ul  121  1.N  7. 
Ill  441  1.SI  1 
221  OM  1.12  0. 
221  IN  l.tt  I 
114  16  1.01  0. 
71  16  1.U  4. 
16  ISS  1.10  4. 
217  06  1.U  4. 

220  60  l.e  4. 
200  SON  2.N  20. 
2U  UM  1.6  4.1 
220  IIN  1.6  M.i 
220  12N  2.6  U.1 
121  16  X.U  0.3 
16  27<  1.6  2.1 


71  4.1  1.15 
17  tt.C  s.n 
W  0.1  2.77 
90  1.7  4.11  , 

IS  w.i  u.n  . 

I  U.7  I.1S  , 
«  4.1  4.K  . 
7  1.1  $.11  . 
1  U.t  14.10  . 
1  1.1  4.SC  . 
9  4.1  S.6  j 
9  I7U  1.$4  j 
9  U.7  tl.M  .3 
3  10.2  S.6  .3 
:  U.0  t.e  .3 
U.7  7.M  .1 
0.7  2.14  .1 
U  2.6  .1 
11.1  U.6  .0 
U.S  1.12  .4 
U.0  M.e  .4 
U.S  f.SI  .« 
11.1  ii.e  .K 
11.f  U.20  .K 
U.S  0.6  .K 
11.0  0.M  .6 
U.S  0.74  .N 
14.0  1.02  .N 
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SAMPLE  SUPERCALC  OUTPUT 


Slope  computations  for  Fisk 


1 

2  DATE 

3  Ifl/lf/77 

<  n/15 

5  J2/H 

6  2/27/78 

7  3/22 

8  4/18 

9  5/17 

10  6/22 

11  7/26 

12  8/31 

13  9/20 

14  10/19 

15  Il/ZI 

16 

17 

18 

19  3/22 

20  4/20 

21  5/24 

22  6/21 

23  7/25 

24  8/22 

25  9/19 

26  10/25 

27  11/20  3.08 

28  12/19  10.76 

29  1/25/80  4.36 


30  2/19 

12.48 

31 

3/20 

•  8.36 

32 

4/23 

11.02 

33 

5/21 

7.59 

34 

6/18 

1.1 

35 

7/50 

.94 

36 

8/21 

.16 

37 

9/17 

.26 

38 

10/21 

.28 

39 

11/19 

2.07 

40 

12/17 

10.07 

41 

1/22/81 

.6 

42 

2/18 

4.4 

43 

3/25 

3.92 

44 

4/22 

1 

45 

5/21 

13.56 

46 

6/17 

13,07 

47 

7/23 

11.96 

48 

8/19 

1.84 

49 

9/24 

9.82 

SO 

10/21 

.59 

51 

11/17 

1.57 

52 

12/16 

5.22 

53 

1/20/82 

2.44 

54 

2/24 

20.24 

55 

5/17 

12.91 

56 

4/20 

7.02 

57 

5/19 

.95 

58 

6/22 

10.12 

59 

7/21 

.92 

60 

8/25 

.13 

61 

9/22 

12.05 

IJ  E  li  f  •;  6  ; 

GLENKONVIUE,  MO 
CASE  HT  287.74S10PE 
2.35  .0001706 

.0001912 
.0001908 
.0001756 
.0003339 
.0001906 
.0001726 
.0001525 
.0001764 
.0001196 
.0001633 
.0001599 
.0001773 
.0002105 
.0002068 
.0001130 
.0002407 
.0002421 
.0002264 
.0001661 
.0001519 
.0001579 
.0001569 
.0001559 
.0001722 
.0002140 
.0001796 
.0002119 
.0001964 
.0002112 
.0001970 
.0001622 
.0001602 
.0001567 
.0001576 
.0001585 
.0001661 
.0002124 
.0001572 
.0001779 
.0001767 
.0001547 
.0002231 
.0002192 
.0002413 
.0001620 
.0002110 
1.28  .0001581 

1.11  .0001676 

1.86  .0001751 

2.74  .0001613 

11-06  .0002401 

6.93  .0002135 

3-82  .0001904 

1.88  .0001561 

4.89  .0002073 

1.49  .0001591 

4.88  .0001244 

5.73  .0002163 


fISi;,  MO  ZERO: 


12/19 

1/25/79 

2/22 


SAGE  KT 
3,17 
7.72 
14.62 
5.55 
20.48 
8.57 
3.84 
.59 
4.31 
4.61 
1.96 
1.3 
4.54 
10.72 
11.58 
MOTHADE 
19.78 
22.36 
16.46 
2.05 
.7 
.86 
.39 


307.460ATE 

10/20/77 


2/28/78 
NOT  MADE 

5/18 

6/21 

8/30 

9/21 

10/20 

11/22 

12/20 

2/21 

3/21 

4/19 

5/23 

6/20 

7/26 

8/21 

9/18 


4.42 

11.37 

4.13 

5.34 

2.78 

1.95 

2,8 

9.93 

2.04 

1.77 

2.91 
5.1 
6.4 

6.12 

10.52 

12.94 

8.92 

1.78 
2.13 
1.57 
1.22 


10/24 

.95 

2.07 

4.72 

1/24/80 

2.46 

2/20 

6.69 

4.44 

5.32 

3.59 

1.29 

1.37 

8/20 

1.02 

1.01 

10/20 

.92 

1.S5 

4.22 

1.4 

2.7 

3/24 

2.37 

4/21 

2.1 

6.42 

6.4 

7/22 

2.63 

8/18 

2.06 

9/23 

4.14 
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c 


:  A 

It  »  II 

CUD* 

1  £  ; 

«  10/19 

2.44 

KOT  MADE 

«  11/15 

2.33 

11/16 

1,76 

64  11/29 

6.22 

11/30 

4.1 

65  1/10/83 

22.07 

1/11/83 

13.72 

66 

i  2/14 

15.24 

2/15 

7 

61 

'  3/15 

7.78 

4.09 

68 

!  4/13 

11.06 

S.41 

69 

'  5/9 

22.21 

5/10 

13.78 

70 

6/14 

6.64 

6/13 

4.34 

71 

7/12 

1.19 

KOT  HADE 

72 

8/8 

.35 

KOT  HADE 

73 

9/13 

-.63 

MOT  MADE 

74 

10/20 

1.96 

KOT  MADE 

75 

11/9 

9.71 

3.65 

76 

1/16/84 

16.84 

8.68 

77 

2/13 

5.32 

9.78 

78 

3/12 

11.52 

5.64 

79 

4/16 

14.22 

7.55 

80 

5/14 

13.52 

7.12 

81 

6/4 

5.76 

3.61 

82 

6/18 

1.11 

1.82 

83 

7/16 

1.45 

KOT  HADE 

84 

8/13 

-.5 

KOI  KAOE 

8S 

9/11 

-.86 

KOT  HADE 

86 

10/22 

12.82 

MOT  HADE 

87 

11/27 

15.98 

13.88 

88 

12/18 

15.09 

12/17 

8 

89 

1/14/85 

22.73 

14.22 

90 

2/11. 

6.7 

6.56 

91 

3/11 

16.73 

10.24 

92 

4/8 

21.66 

13 

93 

5/13 

12.76 

7.09 

94 

6/11 

13.7 

6.92 

95 

7/15 

13.75 

KOT  HADE 

96 

8/12 

6.65 

KOT  KAOE 

97 

9/16 

3.21 

KOT  HADE 

98 

10/15 

.77 

KOT  HADE 

99 

11/18 

14.56 

6.92 

100 

12/16 

15.94 

e.os 

101 

1/21/86 

10.7 

5.77 

102 

2/18 

8.32 

5.36 

103 

3/17 

11.31 

6.27 

104 

4/21 

8.73 

7.63 

105 

5/19 

2.09 

2.91 

106 

6/16 

12.22 

6.4 

107 

7/21 

1.13 

KOT  HADE 

108 

8/26 

-.49 

KOT  HADE 

109 

9/16 

-.17 

KOT  HADE 

110 

10/21 

.63 

KOT  HADE 

111 

11/17 

-.93 

1.12 

112 

12/15 

6.47 

3.68 

113 

1/20/87 

3.17 

2.48 

114 

2/23 

6.8 

3.91 

US 

3/16 

8.05 

4.5 

116 

4/13 

3.21 

2.62 

117 

5/18 

-.54 

1.83 

118 

6/15 

-.8 

1.28 

119 

7/13 

4.28 

KOT  HADE 

120 

8/16 

1.33 

KOT  HADE 

121 

9/15 

.66 

MOT  HADE 

122 

10/13 

-l.OI 

KOT  HADE 

.OOOU85 

.0401614 

.0002332 

.0002323 

.0001945 

.0002107 

.0002338 

.0001829 


.0002141 

.0002316 

.0001268 

.0002127 

.0002192 

.0002170 

.0001817 

.0001579 

.0001759 

.0001597 

.0001567 

.0002703 

.0001813 

.0002227 

.0002345 

.0001650 

.0002177 

.0002357 

.0002109 

.0002201 


.0002273 

.0002293 

.0002048 

.0001884 

.0002057 

.0001729 

.0001570 

.0002122 


.0001468 

.0001870 

.0001695 

.0001878 

.0001933 

.0001687 

.0001441 

.0001465 


179 


i  A  i 

8  ;; 

c  1}  0  ;; 

E  a  f 

i!  c  ; 

123  11/17 

-.8 

2.65 

.0001352 

124  12/7 

■  5.12 

4.1 

.0001723 

125  1/20/88 

14.47 

11.47 

.0001887 

126  2/16 

8.72 

5.52 

.0001904 

127  3/14 

9.21 

5.52 

.0001945 

128  4/11 

11.11 

6.42 

.0002028 

129  5/24 

1.01 

2.3 

.0001531 

130  6/20 

.15 

1.54 

.0001523 

131  7/19 

-.32 

KOr  HADE 

132  8/15 

-1.18 

HOT  HADE 

133  9/12 

-.28 

NOT  HADE 
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SAMPLE  SAM  INPUT  AND  OUTPUT 
Wilhelmina,  November  15,  1977 
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* 

*  HYDRAULIC  DESIGN  PACKAGE  FOR  FLOOD  CONTROL  *  • 

*  CHANNELS  (SAM)  * 

*  + 

*  SEDIMENT  TRANSPORT  CALCULATIONS  * 

*  * 

*  VERSION  3.05  31  March  1993  * 

*  * 

*  A  Product  of  the  .  * 

*  Flood  Control  Channels  Research  Program  * 

*Hydraulics  Laboratory,  Waterways  Experiment  Station* 


Msg  1:  SED.  READING  INPUT  DATA  FROM  FILE  [  sed.in  ]  THIS  DIRECTORY. 


TABLE  1.  LIST  INPUT  DATA. 

TIWILCO  77/11/15 
TF  TOFFALETI.  YES 
TF  YANG.  YES 

TF  ACKER-WHITE.  YES 

TF  COLBY  YES 

TF  TOFFALETI-SCHOKLITSC  YES 
TFMPM(1948).  YES 

TF  BROWNLIE,D50  YES 

TF  TOFFALETI-MPM  YES 
TF  LAURSEN(MADDEN),  1 985  YES 
TF  LAURSEN(COPELAND)  YES 
TF  YANG,D50  YES 

TF  ACKER-WHITE, D50  YES 
TF  MPM(1948),D50  YES 

F#  45678  2345678  2345678  2345678  2345678  2345678  2345678  2345678  2345678 

2345678 

TR  -1 

VE  2.06 

DE  4.64 

WI 154. 

QW  1470. 

ESlllE-5 

WT51.8 

PF  0.5  0.5  100.  .25  86.0  .125  3.4 

PFC.0625  2.0 

$JOB 
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1 

TABLE  5.0  SUMMARY  TABLE:  TOTAL  BED-MATERIAL  SEDIMENT  ' 
DISCHARGE,  TONS/DAY 

Q  I  WATER  I  TRANSPORT  FUNCTIONS 

NO  I  DISCHARGE  1  TOFFALETI.  |  YANG.  1  ACKER- WHITEj  COLBY  | 

I  (HEC-6)  I  (HEC-6)  |  (HEC-6)  | 

1|  1470.1  1455.1  3861.1  22002.|  858 1 


Q  I  WATER  I  TRANSPORT  FUNCTIONS 

NO  I  DISCHARGE  I TOFFALITI-  |  MPM(1948)|  BROWNLIE,  |  TOFFALETI- 1 
SCHOKLITSCHI  (HEC-6)  |  D50  |  MPM  | 

1 1  I470.I  2968.1  187.1  1474.|  1642.| 

Q|  WATER  I  TRANSPORT  FUNCTIONS 
NO  I  DISCHARGE  I  LAURSEN  |  LAURSEN  |  YANG,D50  |  ACKER- 

I 

(MADDEN),  85|  (COPELAND)  |  |  WHITE,  D50 

I 

1 1  1470.1  6326.1  1629.1  3819.| 

5153.1 


Q|  WATER  I  TRANSPORT  FUNCTIONS 
NO  f  DISCHARGE  |  MPM(1948),  | 

D50  I 

1  I  1470.1  194.1 

End  of  Job  PRINTOUT  SAVED  IN  FILE  sed.out 

1 
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SAMPLE  STATGRAPHICS  INPUT 

Observed  water  and  sediment  discharge,  computed  bed 
material  transport,  water  year  1978,  Fisk 


FILE: 

Row 

FISKl 

FLOW 

05/30/96 

OBSQS 

09:58:05  AM 

SANDQS 

TOFF 

1 

303. 

33.54 

7.38 

1 . 

2 

925. 

177,32 

77 . 12 

8. 

3 

2860. 

1436.29 

'  1133.23 

79. 

4 

654. 

63.57 

27.34 

2. 

5 

5600. 

1723.68 

1172.10 

1286. 

6 

1250. 

600.75 

474.59 

14  . 

7 

352. 

233.80 

147.29 

1. 

8 

38. 

7.82 

3.44 

1. 

9 

433. 

122.76 

45.42 

1. 

10 

475. 

374.49 

112.35 

1. 

11 

195. 

70.03 

30.11 

1. 

FILE: 

Row 

FISKl 

YANGH6 

05/30/96 

AWH6 

09:58:05  AM 

COLBY6 

TFSCH 

1 

7. 

l.lOOOOOOOEl 

1. 

1 . 

2 

92. 

1.16900000E3 

1. 

39. 

3 

350. 

1.12280000E4 

130. 

190. 

4 

24. 

4.10000000E1 

96. 

3  . 

5 

3853. 

3.60742000E5 

391. 

2304  . 

6 

83. 

1.65000000E2 

1. 

31  . 

7 

10. 

1.40000000E1 

1. 

1 . 

8 

1. 

.  l.OOOOOOOOEO 

1. 

1  . 

9 

13. 

9.00000000EO 

1. 

2. 

10 

4. 

9.00000000E0 

1. 

1 . 

11 

2. 

l.OOOOOOOOEO 

1. 

1 , 

FILE; 

FISKl 

05/30/96 

09:58:05  AM 

Row 

MPMH6 

BROWN 

TOFMPM 

LRSMAD 

.  1 

1. 

2. 

1. 

1. 

2 

7. 

18. 

15. 

1. 

,3 

34. 

186. 

113. 

1057. 

4 

3. 

8. 

5. 

17. 

5 

133. 

1350. 

1419. 

15170. 

6 

9. 

37. 

24. 

82. 

7 

1. 

2. 

1. 

1. 

8 

1. 

1. 

1. 

1 . 

9 

1 . 

5. 

3. 

1 . 

10 

1 . 

2. 

1. 

1 . 

11 

1. 

1. 

1.  . 

1  - 

FILE:  FISKl 

Row  LRSCOP 

05/30/96 

YANG50 

09:58:05  AM 

AW50 

MPM 

1  4 . 

7. 

4. 

0. 

2  82. 

69. 

123. 

8. 

3  245. 

318. 

393. 

41. 

4  18. 

22. 

12. 

2. 

5  2246. 

3116. 

13580. 

145. 

6  38. 

78. 

46. 

11. 

7  5. 

10. 

5. 

1. 

8  1. 

1. 

1. 

0. 

9  5. 

15. 

7. 

1 . 

10  3. 

18. 

78. 

2. 

11  3. 

2. 

1. 

0. 

Page 


Page  2 


Page  3 


Page  4 
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SAMPLE  QUATTRO-PRO  INPUT 
Partial  entropy  calculation  file,  Fisk,  MO 
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INTERVAL 

YANG  HEC6 
FISKARTH 

101 

OBSV 

100  FREQUENCY 

PROBABILITY 

ENTROPY 

1 

54 

0.53465346535  0.33476603361 

2 

7 

0.06930693069  0.18499477797 

3 

5 

0.0495049505 

0.14879616854 

4 

6 

0.05940594059 

0.16772441867 

5 

7 

0.06930693069 

0.18499477797 

6 

4 

0.0396039604 

0.1278743032 

7 

3 

0.0297029703 

0.10445073945 

8 

6 

0.05940594059 

0.16772441867 

9 

1 

0.0099009901 

0.04569426254 

10 

1 

0.0099009901 

0.04569426254 

11 

1 

0.0099009901 

0.04569426254 

12 

0 

0 

0 

13 

0 

0 

0 

14 

1 

0.0099009901 

0.04569426254 

15 

0 

0 

0 

16 

0 

0 

0 

17 

0 

0 

0 

18 

1 

0.0099009901 

0.04569426254 

19 

0 

0 

0 

20 

0 

0 

0 

21 

1 

0.0099009901 

0.04569426254 

22 

0 

0 

0 

23 

0 

0 

0 

24 

0 

0 

0 

25 

1 

0.0099009901 

0.04569426254 

26 

0 

0 

0 

27 

1 

0.0099009901 

0.04569426254 

28 

0 

0 

0 

29 

0 

0 

0 

30 

0 

0 

0 

31 

0 

0 

0 

32 

0 

0 

0 

33 

0 

0 

0 

34 

0 

0 

0 

35 

0 

0 

0 

36 

0 

0 

0 

37 

0 

0 

0 

38 

0 

0 

0 

39 

1 

0.0099009901 

0.04569426254 

40 

0 

0 

0 

101 

1 

1.83257400095 

FREQUENCY 

50 

7 

5 

4 

7 


FISK10  93 

Y  PROBABILITY  ENTROPY 
50  0.5376344086  0.33364327297 
7  0.0752688172  0.19469704741 

5  0.05376344086  0.15715922477 
4  0.04301075269  0.13532495192 
7  0.0752688172  0.19469704741 
3  0.03225806452  0.11077378079 
2  0.02150537634  0.08256886694 

6  0.06451612903  0.17682838864 
1  0.01075268817  0.04873762896 
1  0.01075268817  0.04873762896 
1  0.01075268817  0.04873762896 
0  0  0 
0  0  0 
1  0.01075268817  0.04873762896 


1  0.01075268817  0.04873762896 


1  0.01075268817  0.04873762896 


1  0.01075268817  0.04873762896 
0  0  0 
1  0.01075268817  0.04873762896 


1  0.01075268817  0.04873762896 
0  0  0 
33  1  1.82433124146 
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herein  will  define  the  amount  of  an  existing  sediment  data  set  that  provides  the  most  information  and  also  evaluate 
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c.  A  procedure  has  been  developed  that  will  allow  evaluation  of  the  adequacy  of  existing  data  sets  for  selecting 
an  appropriate  sediment  transport  equation. 

d.  A  set  of  sediment  data  samples  has  been  assembled  that  is  unprecedented  in  scope  and  completeness,  and 
available  in  digital  format. 

€.  In  sight  has  been  provided  that  will  answer  the  question,  “How  much  data  are  needed  for  a  riverine  sediment 
study?” 


